


Institutional Archive of the Naval Postgraduate 


schoal 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


1966 


1. Thesis and Dissertation Collection, all items 


Non-equilibrium condensation of carbon 


dioxide in superso 


Duff, Karl M. 


http://ndl.handle.net/10945/9591 


nic nozzles. 


Downloaded from NPS Archive: Calhoun 


DUDLEY 
KNOX 
LIBRARY 





http://www.nps.edu/library 


Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarly author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 








PYG EY RNC b icghe 
NAVAL POSTGRADUATE. stpau% 
MOMTFREY CA $15 44654)"' 


HON~BQUILEBRIUM CONDENSATION OF CARBON DIOXIDE 
IN SUPERSONIC NOZZLES 


by 
Lt. Karl a en USN 


5S. 3., Massachusetts Institute of Technology 
(1957) 


S. M., Massechusetts Institue of Technolosy 
(1964) 


Cambridge, Massachusetts 


Subuitted in Partiel Fulfillment of 
the Requirements for the Desree of 


DOCTOR OF SCIENCE 
in 
MECHANICAL EHGINEFRING 


at the 
MASSACHUSETTS INSTITUTE OF TECHROLOGY 


January 1066 





HOW-KQUILIBARING CONDENSATION OF 
TARSON DIOKIDS Th SUPERBORIC NOZZLES 


by 


Karl “%. baff 


Submitted to the Repartment of Mechanical Ingineering on January 16, 
1966, in partiel fulfillment of the requirements for the Degree of Doctor 
of Selence. 


ABSTRACT? 


Yests involving the concensation of C0, in three supersonic nozsles 
have been conducted. Late of streanwise pressure variation and local 
eonditions at anset of condensation heve been obteined and show that 
rate of expansion, ag determined vy nozsle geometry, leas a significent 
influence on the supersaturation obtained prior to sondensetion. 


Interferometric density measurements for non-condensing flow 
verified the applicebility of l-~dimensional vas dynamics to the expen- 
siong, and the need for consideration of departure from perfect raz 
behavior of CO Non-condensing density--pressure meagurenments and 
total pressure reasurements subatantiated the use of the Plank Zquation 
of State in the low temperature region of testing. 


Classical nucleation and drop cng theory, as apclied in this 
study tc nozzle flow, will predict the conditions of suyersaturetion 
at onset of condensation reasonably well , provided amall corrections 
to the flat film liquid surface tension of CO. are applied. These 
corrections amount to increases of from 177 to 328%, sepending uroen 
Values of thermal and masse seccommodation coefficients whieh are also 
applied. 


The applied theory does not suceassfully predict the experimental 


Streamvise veriation in pressure caused ty condensation, mal the extent 
of effeeta of nozzle ceometry noted experimentally These falladngs are 


thought to indicate inadequacies in drop growtha theery, or in its 
application. 
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NOMENCLATURE 


A Nozzle Flow Area 

os Specific Heat at Constant Pressure 

C. Specific Heat at Constant Volume 

gz Number of Molecules in a Droplet Cluster 

aX Number of Molecules in a "Critical" Cluster 
h Specific Enthalpy 

Se Specific Enthalpy of Vaporization 

H Total Enthalpy 

J Number of Nuclei Produced per Unit Time per Unit Volume 
k Boltzman Constant 

m Molecular Weight 

m' "Effective'’ Molecular Weight 

m Mass Flux, Particles/Unit Area-sec. 

M Mach Number 

n Number of Vapor Molecules/Unit Volume 

n Index of Refraction of a Vapor in State i 


N(g) Equilibrium Size Density Distribution of g-molecule Clusters 
PS op Static Pressure 

i Stagnation Pressure 

ae Stagnation Pressure behind a Normal Shock 


Pov D Flat Film Fquilibrium Vapor Pressure 


Py Droplet Equilibrium Vapor Pressure 
Q Heat 

va Droplet Radius 

r* Helmholtz "Critical" Radius 

R Gas Constant 


ae MeLecetve. Gas Constant 
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hntropy 

Fringe Shift, Number of Fringes 

Time 

Static Temperature 

Stagnation Temperature 

Drop Temperature 

‘Temperature of Molecules Reflected from a Droplet 
Streamwise Velocity 

Internal Energy of a Droplet 

Specific Volume 

Total Volume 

Streamwise Variation of Distance in Nozzle 
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Boundary Layer Thickness 
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Quasi-Equilibrium Droplet Density Distribution 
Angle, Radians 

ee 

Wave Length of Interferometer Light Source - 5461 a° 
Mass Fraction of Condensed Moist.ure 
Condensation (Mass Accommodation) Coefficient 
Local Density 


tagnation Density 





vil 
Droplet (Cluster) Density 
Drop Surtace Tension 
Liquid Flat Film Surface Tension 
aby AL 
Free Energy per Molecule (Subscripts; L = Liquid, V = Vapor) 


Total Free Energy of g Molecules 
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1 INTRODUCTION 
In recent years, increasing attention has been given to the problem of 
homogeneous nucleation and condensation of supersaturated vapors. 
Particular interest has arisen in connection with high temperature, 
low vapor pressure fluids, capable of high Carnot efficiencies in power 
plant application. The possible space application of systems employing 
alkali metals, and recent reports of experiments involving the flow of 


43,54 


sodium and potassium through nozzles underlines the need to understand 
the condensation behavior of these gases well enough to provide a basis for 
optimum turbine and cycle design. 

Early investigations in cloud chambers showed that, after initial 
expansions and settling had removed dust particles and other nuclei, sub- 
sequent expansions into the saturated region would result in the spontaneous 
formation of a condensation cloud at some limiting expansion ratio. In the 
absence of other nucleation sites, it was concluded that such clouds were 
the result of homogeneous nucleation of the super-saturated vapor. Similar 
condensation occurs in vapors expanded in nozzles, where rates of expansion 
are generally much more rapid than those of cloud chambers. 

Despite extensive literature on the subject, there exists little agree- 
ment as to the correct formulation of nucleation and drop growth rates, and 
within the theoretical framework there exist many other uncertainties regard- 
ing the fluid properties pertinent to numerical evaluation. Experimental 
work is relatively limited in terms of range of test conditions, and numbers 
of fluids tested. Although a variety of gases have been studied in cloud 
chamber expansions, applicable and informative studies of nozzle expansions 
are somewhat sparse outside of steam and airborne water vapor for which both 
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cloud chamber and nozzle data is extensive. This data as well 


as the cloud chamber data for several other vapors tends to substantiate 
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formation of a condensation cloud at some limiting expansion ratio. In the 
absence of other nucleation sites, it was concluded that such clouds were 
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condensation occurs in vapors expanded in nozzles, where rates of expansion 
are generally much more rapid than those of cloud chambers. 

Despite extensive literature on the subject, there exists little agree- 
ment as to the correct formulation of nucleation and drop growth rates, and 
within the theoretical framework there exist many other uncertainties regard- 
ing the fluid properties pertinent to numerical evaluation. Experimental 
work is relatively limited in terms of range of test conditions, and numbers 
of fluids tested. Although a variety of gases have been studied in cloud 
chamber expansions, applicable and informative studies of nozzle expansions 
are somewhat sparse outside of steam and airborne water vapor for which both 


7,62,63,66 & 67 


cloud chamber and nozzle data is extensive. This data as well 


as the cloud chamber data for several other vapors tends to substantiate 
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so-called "classical nucleation theory, as opposed to revised models 
of nucleation theory offered by more recent authors, as a better means of 
predicting condensation behavior. However, there is doubt regarding the 
extent to which cloud chamber data may be applied to the determination of 
nucleation rates, and too little data outside of water vapor has been obtained 
for flow in nozzles to allow a comprehensive test of theory over a variety 
of gases. 

The purpose of this work has been to extend the limits of existing 
experimental data involving condensation in nozzles, through experimentation 


with CO,, to study by means of a computational program the predictions of 


a? 
existing theory and make some judgments regarding the importance of the 
various physical parameters and uncertainties pertaining thereto, and to 


reduce these uncertainties to a point where a valid comparison of existing 


theory and experimental results can be made for CO, « 





II, NOZZLE FLOW PHENOMENA 

Converging-diverging nozzles lend themselves particularly nicely to 
the study of condensation behavior of gases. Following are some of the 
aspects of gas behavior in nozzle flow and features which recommend it as a 
means of experimental study of condensation. 

A. Isentropic Effects 

It is general knowledge that for sufficiently high pressure ratios, an 
ideal gas will expand isentropically in a converging-diverging nozzle, 
becoming sonic at the point of minimum flow area and supersonic downstream, 
with continuously increasing velocities and decreasing pressures, densities 
and temperatures all along the expansion. A non-ideal gas behaves similarly, 
and depending upon the stagnation conditions, the extent of the expansion, 
and the properties of the vapor being tested, the isentrope may approach 
saturation conditions of pressure and temperature. As illustrated on the 
T-S diagram of Figure 1, it is possible for conditions to exist for certain 
gases when the opposite may be true and an isentropic expansion will result 
only in increasing the degree of superheat of the vapor (i.e. as with many of 
the FREONS). However, in most cases it is a simple matter to devise stagna- 
tion conditions which result in the rapid arrival of the gas to conditions 
of supersaturation during its isentropic expansion. For a supersonic nozzle 
a few inches in length expansion time is of the order 107" seconds. 

Be Wall and Boundary Layer Effects 

In all real gases there are departures from ideal non-viscous behavior 
which must be considered. Fortunately, in nozzle flow these effects are 
easily accounted for, provided fully developed flow does not occur. Viscous 
effects are accounted for by considering the nozzle flow to be an isentropic 
expansion contained between wall boundary layers. The boundary layers are 


greatly influenced by the shape of the pressure gradient and are generally 
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depressed to a very small mimimum value near the most critical point, the 
geometric throat. Further, viscous and recovery effects on the nozzle walls 
will keep the wall temperature close to stagnation value in the absence of 
very strong wall cooling or heating, and, at decreased downstream pressures, 
is adequate guarantee that there will be no condensation of liquid on the 
walls which could interfere with the mid=stream process. 

In experimentation where nozzle size might be limited by mass flow or 
other restrictions, considerations of Reynold’s Numbers and boundary layer 
thicknesses which may be expected in the throat region are important since 
large boundary layers in relation to nozzle size could invalidate the 
assumption of an isentropic region of flow. A means of estimating the 
boundary layer thickness in the throat region of a supersonic nozzle was 
devised (Appendix E) and showed good agreement with interferometric photos 


that the boundary layers involved in the CO, experiments of this investiga= 


2 
tion were about .006" at the nozzle throat. In all nozzles tested this was 
less than 11% of the throat width. Additional discussion of boundary layer 
treatment in this investigation is included later. 

The isentropic region between boundary layers may be described by one= 
dimensional gas dynamics, and, by means of measurement of some vapor 
property (such as static pressure) along the nozzle, a complete time- 
history of the expansion may be obtained. 

Condensation in a nozzle expansion has been shown by many investiga- 
tions to have a large influence on the thermodynamic properties of the 
fluid. In particular, the static pressure records large variation from the 
non=condensing profile as a result of the heat released by moisture forma- 
tion. For the addition of heat to an ideal gas, it may be shown that 


dP _ yM¢ (aA dQ_j 


Pp Toe “Ae te cor 
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For water vapor, it is easily estimated that at M = 2, a formation of 
1% moisture at constant area is accompanied by about 7% variation in static 
pressure. As is later shown, it is possible to detect the onset of condensa- 
tion at about .1% presence of moisture through simple static pressure 
measurements. 

It may be seen from the above equations that pressure variation is 
increased by establishing conditions such that condensation occurs with 
Mach Number close to unity. 

Willmarth and Nacagecenaes and Arthur,> have shown that the ratio E/ES 
plotted against the streamwise variation of ae is an even more sensitive 
indicator of condensation, where P is the measured static pressure, Pe is 
the measured total pressure behind a normal shock, and Po is the original 
stagnation pressure. 

An important consideration in any system designed to test homogeneous 
nucleation theory is the possible effects of contamination. As will be shown 
later, nozzle expansions from superheated to supersaturated state proceed in 
such an extremely short time that spontaneous nucleation will produce nuclei 
many orders of magnitude greater in number than those which could be present 
in the form of dust or other particles and the amount of moisture formed on 
such contaminants is entirely negligible. As will also be mentioned later, 
an exception to this is the case of a pre-condensing vapor. 

The proceeding observations indicate that, with fairly simple instru- 
mentation of a nozzle system, it is possible to obtain a complete history of 
@ gas expansion and any condensation occurring within, with the reasonable 
assumption that no outside influences exist which would invalidate testing 
of homogeneous nucleation and drop growth theory in the same system. 


67 Test Vapor Considerations 


The selection of a vapor for testing is subject to the following 
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considerations. First, a large number of vapors are toxic, reactive to 
instrumentation materials, or do not exist as vapors at conditions most 
convenient for laboratory testing. Second, the state relations and other 
properties of many are so poorly defined, that application to any meaningful 
test of theory would be ill-founded. Neither of these considerations would 
be restrictive if the condensation data for a particular vapor were parti- 
cularly needed, or if a test of theory were incidental to the experimental 
data. Third, two particular properties are desirable in a test vapor, both 
of which are more common to gases of low molecular weight. The first has 
been mentioned previously, and involves the shape of the coexistence line in 
relation to the isentrope. A preferred test gas has a coexistence line 
having rapid convergence with the isentrope. Also, the measure of heat 
release resulting from condensation compared with that needed to alter its 
vapor properties is important. This may be expressed through the parameter 
hep /C,» for which typical values for various gases are presented in Table l. 
Variations in the gas properties caused by condensation are increased with 


larger values of this parameter. A more applicable use of these values is, 


as seen from Equation 1, the form h 


fe oe at the test conditions of interest. 


For vapors such as alkali metals, higher values of ae are decreased in 
Significance due to the high temperatures of the gas. 

CO, was selected for this study because it satisfies all the foregoing 
criteria and in addition may be obtained in high purity at relative low 
cost, thus allowing its use in a system not necessitating recovery after 


expansion. Also, such physical properties as CO, vapor pressure, liquid 


2 
surface tension, and solid and liquid density, necessary in the application 


of theory, are well established and readily available in numerous publica- 


tions of physical property data. 





Discussion of the application of these properties, and details of the 
theoretical mechanism and behavior of condensation are presented later. 
The question could be raised concerning possible non-equilibrium of 


energy modes of CO, during such a rapid expansion. Kenereat tae has shown 


a 
that this would not be a problem with CO, in this case since 1)vibrational 
energy modes are not activated at the low temperatures of experimental test 
conditions and 2)the relaxation time constant of vibration is considerably 
less (aionccaes than that of the gas expansion. Conacwe has recently 
formulated the pressure=-temperature dependence of CO, vibrational relaxa- 


tion time and his expression also yields times of the order 107° sec. in 


the region of experimental testing. 





TII, EXPERIMENTAL INVESTIGATIONS 


ee Test Apparatus 


Industrial Welding CO, having very high (moisture content less than 


2 
.0045% by weight) purity was introduced via a manifold from a bank of nine 
50 1b cylinders to a needle valve used for regulation of flow rate (See 
Figure 2), Thence the flow entered a stagnation tank and then the test 
nozzle. Exhaust was to atmosphere in earliertests and later to a steam 
ejector system maintaining apour 23" Hg. of vacuum. This modification 
facilitated lower minimum stagnation pressures without the interference of 
compression shock systems in the supersonic test section. Nine cylinders 
were employed in order to obtain a large mass flow source providing nearly 
constant conditions to the nozzles for extended periods of flow, despite 
the quasi-steady nature of the "blow-down" method. 

The stagnation tank contained a calibrated pressure gauge and, extend- 
ing into the flow just upstream of the nozzle entrance, a copper-constantan 
thermocouple using an ice reference. A 100 inch 25=tube mercury filled 
manometer board was used in connection with static wall taps to cetermine 
pressure readings along the nozzle. Pressure data was recorded by 
polaroid photograph of the manometer board. Test time durations were of the 
order of 1 minute. 

The stagnation tank had a removable end plate and contained screens 
which allowed the packing of large amounts of steel wool inside a nylon 
filter bag upstream of the nozzle entrance. (Figure 5) When high stagnation 
temperatures were desired (i.e. non-condensing tests) an electric heating 
tape wrapped around the stagnation tank was energized and the steel wool 
used for heat storage and thermal inertia to decrease the rate of the 
temperature transient. When low temperatures were desired, the steel wool 


was removed and the heating tape not used. 





Bs Test Nozzles 

Two test nozzles were employed for experimental measurements, 
henceforth referred to as Nozzle I and Nozzle II. 

Nozzle I appears schematically in Figure 7 and in photos in Figures 6 
and 9. It featured two removable parallel side walls which were plates 
of optical glass, used for interferometric measurements. By exchanging 
one of the plate; for an identical metal one containing pressure taps, 
static pressure profiles could be obtained and the flow observed simulta- 
neously through the remaining glass wall. Slotted "O-rings" provided 
sealing of the side walls. The nozzle was made wide in comparison to its 
height in order to minimize any contribution to error caused by boundary 
layer effects of the side walls on the interferometric measurements. This 
error was later estimated to be less than .2% of the measured densities. 

The nozzle shape was also intended to provide nearly 2-dimensional flow, 
and in fact, interferometric photos to be discussed later show the flow to 
be effectively 1-dimensional between boundary layers. A region of isentropic 
flow was thus obtained between the boundary layers, which were of thickness 
=~.005 = .007 Jab-the throat | fppendi x. fh). 

Pressure Taps in Nozzle 1 were spaced 1.2 “apart. and tater durine the 
program, additional taps, resulting in .1" spacing, were added in the test 
section region. 

Nozzle II is illustrated in Figures 8 and 10 and featured hinged 
upper and lower nozzle walls which could be placed at arbitrary angles 
between the fixed parallel side walls. This facilitated variable expan- 
sion rates in the test section, depending on the divergence angle of the 
test section, and it was hoped to thus establish a “best” angle of 
divergence to obtain the most informative pressure profiles. Since 


pressure variations due to heat addition approach maximum when Mach 
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Number approaches 1.00, it was expected that this result would be obtained 
for a very low angle of divergence. 

Static pressure taps located in the side wall were located .25 inches 
apart for early runs, after which a new test plate with taps .1 inches 
apart in the test region was substituted. 

C. Experimental Procedures 

i Pressure-Density Measurements 

Tests were conducted using two persons and the following procedure. 

1) The needle valve was opened until the approximate value of 
desired stagnation pressure was obtained. 

2) One observer maintained a null-balance on the thermocouple- 
potentiometer circuit, while the other, manning a tripod-mounted polaroid 
camera, observed the manometer board. 

3) When the manometer board reached a steady state, a polaroid 
photo was taken of the manometer board by the second observer while the 
first recorded the stagnation gauge pressure and terminated the null- 
balancing procedure on the thermocouple potentiometer. 

4) The reading of the potentiometer at time of photograph was 
converted to a stagnation temperature reading, and the polaroid display of 
the pressure data developed, and all pertinent data involving the test 
run affixed to it (including local barometric pressure). 

Interferometer runs were made in a Similar manner, except for the use 
of 5" x 7" sheet film in the built-in camera of the interferometer, and 
requiring only the potentiometer null-balance prior to recording data. 

Discussion of interferometry theory and the methods applied in this 


investigation plus related references are presented in Appendix G. 
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2% Stagnation Pressure Measurements 

The stagnation pressure behind a normal shock was measured by means 
of a .030" O.D. pitot tube extending upstream from a micrometer mounted on 
the elbow downstream of the nozzle.(See Figure 6) All measurements were 
made on Nozzle I, by visually aligning (with the aid of lighting, and 
machinist's templates) the tip of the tube with the center of the corres- 
ponding wall tap with which static pressure was measured. Since the 
location of the shock was actually located a small distance upstream, a 
correction was made to the static pressure reading corresponding to this 
distance, (about .63 tube diameters at Mach 1.65 and about .30 tube dia- 
meters at Mach 2.45, decreasing further as Mach Number enero ee 
Upstream stagnation pressures were measured with the calibrated gauge in 
the stagnation tank. As in the case of the static pressure profile measure- 
ments, photographs of the manometer board were taken simultaneously with 
stagnation temperature and pressure measurements, after steady pressure 
readings and a null potentiometer reading were obtained. 

Total pressure measurements were made with all manometer tubes dis- 
connected except the following 1)pitot tube 2)static wall tap at tip of 
pitot tube 3)throat static wall tap (recorded only for supplementary 
purposes). 

Discussion and estimates of possible experimental error are contained 
in Appendix. 

Dy Non-Condensing Measurements 

lie Real Gas Behavior 

It became necessary to consider the departure from perfect gas behavior 

of CO, near the very outset of this investigation, since one of the objects 


of experimental work was to compare the agreement of measured density pro- 


files with measured pressure profiles. Failure to consider the proper 
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aS Stagnation Pressure Measurements 

The stagnation pressure behind a normal shock was measured by means 
of a .030" O.D. pitot tube extending upstream from a micrometer mounted on 
the elbow downstream of the nozzle.(See Figure 6) All measurements were 
made on Nozzle I, by visually aligning (with the aid of lighting, and 
machinist's templates) the tip of the tube with the center of the corres- 
ponding wall tap with which static pressure was measured. Since the 
location of the shock was actually located a small distance upstream, a 
correction was made to the static pressure reading corresponding to this 
distance, (about .63 tube diameters at Mach 1.65 and about .30 tube dia- 
meters at Mach 2.45, decreasing further as Mach Number YS ere ee 
Upstream stagnation pressures were measured with the calibrated gauge in 
the stagnation tank. As in the case of the static pressure profile measure- 
ments, photographs of the manometer board were taken simultaneously with 
stagnation temperature and pressure measurements, after steady pressure 
readings and a null potentiometer reading were obtained. 

Total pressure measurements were made with all manometer tubes dis- 
connected except the following 1)pitot tube 2)static wall tap at tip of 
pitot tube 3)throat static wall tap (recorded only for supplementary 
purposes). 

Discussion and estimates of possible experimental error are contained 
in Appendix I. 

Dp, Non-Condensing Measurements 

ilms Real Gas Behavior 

It became necessary to consider the departure from perfect gas behavior 

Or CO, near the very outset of this investigation, since one of the objects 


of experimental work was to compare the agreement of measured density pro- 


files with measured pressure profiles. Failure to consider the proper 
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p-v-T relation for reduction of interferometric data resulted in consider= 
able error in the density profile. It also introduced doubt as to how the 
proper density ratios should be converted to the corresponding isentropic 
pressure ratios for purposes of making comparison between measurements. 


Existing tabulated properties of CO, also showed that the degree of 


ee 


departure from perfect gas behavior for CO (as reflected by the compressi- 


2 


IZ , : 
bility Tacvor, 2° = aa increased markedly near the saturated region, and 


that the values of compressibility factor could easily be as low as .9 
during an expansion. Also, existing tables did not extend to pressures and 
tefmperatures low enough to make adequate determination of what the real 


gas behavior of CO, is. 


Since a meaningful test of the nucleation and drop growth predictions 
which were to follow must rest on the most precisely accurate information 
available regarding the gas dynamics, and since, as will be shown, the pre- 
dictions are quite sensitive to errors in this respect, it was decided to 


establish the real gas behavior of CO, through use of an applicable equa- 
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Unfortunately, though the literature abounds with work, much of it 


et sole econo. 


recent, on high pressure and temperature properties of CO,» 
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it was possible to find only one equation of state that is 


applicable in the lower P-T range of these experiments, that being 


NT | 


Plank's 1929), which fits the data of Lowry and ieeon and Maass and 


yy 
Mennie : very well. Comparison of Plank's Equation of State and of their 


: : i . : : Py 
data is shown in Figure 1, showing maximum disagreement of about .1% of Rr 


For comparison, the equation of State of Martin & Hou, and Benedict, Webb, 


and Rubin,applied well outside their applicable region to the area of 


experimental interest, vary by as much as 10% from Plank's. 
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Plank's equation is statedand was applied as follows: 
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The right hand term may be evaluated from the equation of state and 
hence it is possible to determine enthalpy at any desired pressure and 
temperature. To determine an isentrope, equation (3) is used and the change 
in h calculated for an arbitrary change in p. Thence, knowing new values 
of p and h, the corresponding value of T may be determined, either through 
the use of pre-calculated enthalpy tables, or by iterative integration of 
equation (4). It is then a simple matter to calculate the new value of 
density from equation (2) and to thus determine a series of pressure, 
temperature, and density ratios along an isentrope for any given set of 
stagnation conditions. 

Table 2 contains values of enthalpies calculated in this manner. 
Agreement within .1° to .3°K occurs at the juncture of those presented 
by Doneee 

In addition to the comparison of pressure and density measurements, it 
was decided to apply an additional experimental test to the ability of the 


Plank Equation of State to predict the behavior of CO, gas dynamics in 


a 


regions of low pressure and temperature yet below its known region of 


applicability, (O0°C to ~70°C). Measurement of the normal shock pressure 
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characteristics were made and compared with the predictions of "real" gas 
dynamics. For this purpose, it was necessary to include a computation 
scheme for solution of the continuity, momentum, energy, and state equa- 
tions across a normal shock and to return along a new isentrope to the new 
stagnation conditions. A listing of this computation procedure is given in 
Appendix H. 

2. Pressure-Density Measurements 

Photographs of two interferometer fringe patterns for non-condensing 
flow are seen in Figures 11 and 12 and the resulting densities given in 
Table 3 and 4. The photographs are marked by an extended region of near 
verticle interference fringes between the upper and lower nozzle walls, 
indicating that this region is in l1-dimensional flow, and that the boundary 
layers are quite small (approximately .006" at the throat). 

Results of the two sets of measured density and the corresponding 
pressure profiles of Nozzle I are shown in Figures 15a, 15b and l6a, 16b. 
Each comparison has been illustrated in two ways. Figures 15a and 16a show 
the density and pressure profiles in relation to nozzle position. Pres- 
sures have been directly measured and plotted at their tap locations. 
Measured density ratios have been converted to their corresponding pressure 
ratios through use of the real gas isentropic values based on the experimen- 
tal stagnation conditions. 

Figures 15b and 16b show the same data plotted in somewhat different 
form, illustrating the theoretical isentropic values of P/P VS. p/d . and 
those values of P/P and p/e which were measured at identical points in 
the nozzle. On Figure 16b, the corresponding perfect gas conversion from 


density to pressure (based on zero pressure value of Cc is also shown. 
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Figures 15 and 16 show very good agreement between measured pressures 
and densities, and those predicted from the Plank Equation of State. The 
fact that the conversion of density to pressure is not very significantly 
affected by the use of the equation of state rather than perfect gas 
relations, is indicated by the close proximity of the perfect gas line in 
Figure 16b, where the applied ratio of specific heats is the zero pressure 
value at stagnation temperature, y = 1.28. However, it should be emphasized 
that even this use of the perfect gas approximations is preceeded by the 
use of real gas p-v-T relations in the reduction of interferometric fringe 
shift data to obtain values of density. The only noticeable departure from 
exact real gas agreement consistent to both tests occurs near the throat 
region where the measured pressures areslightly higher (= 1.8" Hg) than 
those values based on measured density. This departure is concluded to be 
due to departure from 1-dimensional flow over the boundary layer of the 
lateral side walls in this region and an approximate numerical estimate of 
this effect in good agreement with the noted error is made in Appendix F. 

Attempts to apply the interferometer to direct study of condensation 
were frustrated by increasing fluctuation and blurring of the interference 
fringes as stagnation pressure was increased. This effect was apparently 
due to thermal turbulence. Above pressures of 85 psia, complete blurring 
resulted upstream of the throat, although the interference bands again 
became visible downstream, where acceleration of the flow dampens the 
effects of turbulence. Consequently, no density profiles involving conden- 
Sation were obtained although photographs of the condensation ice cloud 
showing the interference fringes within were obtained. 

3. Normal Shock Total Pressure Measurements 
A more sensitive test of the Plank Equation of State predictions 


versus those of perfect gas was made by applying each to the gas dynamics 
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of normal shock behavior and comparing experimental results with theoretical 

predictions based on the experimental stagnation conditions of each test. 
Results of experimental normal shock measurements are shown in 

Figure 17, using the display of Willmarth and Nagamatsu, and Arthur, who 
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made similar measurements for Nitrogen and Helium. The predictions 
computed with Plank's Equation of State are shown in Figure 18. Test 
series No. 2 was taken with a great deal of regard for better accuracy in 
pitot tube positioning and establishment of steady state readings prior to 
taking data, with a resulting decrease in the amount of scatter. One- 
dimensional perfect gas predictions are shown for reference in the lines of 
constant y, the appropriate zero pressure value for most of the test con- 
ditions being about y = 1.28. It is seen that the real gas predictions are 
considerably more accurate than those of perfect gas in predicting the 
normal shock characteristics of CO,» although there still exists some 
departure of the experimental results from those predicted with the Plank 
Equation of State. The experimental results indicate behavior very close 
to that of a perfect gas having a value of y = 1.40. However, this is 


misleading, since treatment of CO, as a perfect gas, using this value of 


2 
ratio of specific heats, does not successfully correlate the pressure- 
density measurements, or correctly predict local temperatures along an 
isentrope. In terms of the measured properties, the experiment indicated 
about 2% larger values of stagnation pressure behind the normal shock than 
those predicted. 

Despite the above mentioned lack of complete agreement between Plank's 
Equation of State and measured normal shock pressure characteristics, it 


"real gas" treatment of the gas dynamics is a con- 


appears that such a 
siderably better description than perfect gas and so warranted application 


to the prediction of local conditions along the isentrope. 
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he. Influence of Equation of State Upon Isentropes 


“he influence of real gas treatment to the computation of CO, isentropes 


2 
has quite a significant effect upon the values of local temperature calcu- 
lated. Figures 19 and 20 compare two typical perfect gas isentropes with 
those predicted using Plank's Equation of State. These curves indicate that 
an error of 5°K could easily occur by attempting to describe the isentrope 
by a constant ratio of er heats, and that a decrease in stagnation 
tempereture could cause errors of possibly 10°K if one isentrope were applied 
to all experimental test conditions. 

The resulting uncertainty in predicted nucleation rate, depending upon 
location within the supersaturated region, can be seen by referring to tne 


classical nucleation rate profiles for CO, illustrated in Figure 21. 


2 


Errors of 15 orders of magnitude could easily occur. 


Ba Condensation Measurements 
ike Pressure Profiles and Local Conditions at Onset of 
Condensation 


Results of experimental measurements involving condensation in 
Nozzle I, pressure tap plates 1 & 2 are shown in Figures 22 and 23. 
Results of Nozzle II, pressure tap plate 1 are shown in Figure 24, and of 
the same nozzle using pressure tap plate 2 in Figures 25 through 30. 
Figures 2h, 25, 26, 27 and 28 are for Nozzle II at an included divergence 
angle of approximately 2.2° while Figures 29 & 30 are for Nozzle II at an 
angle of 0.80°. For convenience and due to manometer board limitations, 
pressure ratios are expressed in terms of the pressure at the tap nearest 
the throat, but may be easily converted to total pressure ratio. Data 
for each Figure is adjoined in Tables 5 through 13. 

All experimental profiles showed a marked occurrence of pressure 


variation from non-condensing values at the point of occurrence of 
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condensation (verified visually in Nozzle I). By noting the pressure ratio 
at which condensation occurs, it is possible to determine the local condi- 
tions at occurrence using the isentropic relations for the equation of 
state previously discussed. Based on the range of applicability of Plank's 
Equation of State and comparison with the tables of Dine 2 it is estimated 
that temperatures in the vicinity of 190°K are accurate to within .3°K and 
temperatures in the region of 160°K accurate to perhaps 1° or 2°K. A plot 
of these points for all three nozzle geometries tested is shown in 
figure 3. 

Several characteristics of the condensation profiles and the local 
conditions at point of occurrence may be noted. 

First, the shape of the profiles contain no sharp minimum at the 
point of occurrence of condensation such as those noted in the condensation 
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of pure water vapor. .6. it is not appropriate to speak of condensa-— 
tion shock in this case, since the formation of moisture occurs along an 
extended region. 

Second it appears that pressure level (or temperature level, or some 
dependent property) may influence the shape of the profile, i.e. the rate 
of formation of moisture, to some extent since the shapes of the condensa- 
tion profiles appear to be more "Spread out" for high stagnation pressures 
than for low stagnation pressures. See Figure 27. These effects are not 
explainable in terms of gas dynamics influences, and the cause is not 
apparent. 

Third, a very definite influence of nozzle geometry upon the point of 
onset of condensation occurs, as seen in Figure 31. The rapid rate of 
nozzle divergence and expansion of Nozzle I results in higher supersatura- 
tion ratios and relatively greater delay of condensation. Progressively 


decreasing the rate of expansion through decrease of nozzle divergence 
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condensation (verified visually in Nozzle I). By noting the pressure ratio 
at which condensation occurs, it 1S possible to determine the local condi- 
tions at occurrence using the isentropic relations for the equation of 
state previously discussed. Based on the range of applicability of Plank's 
Equation of State and comparison with the tables of Dee it is estimated 
that temperatures in the vicinity of 190°K are accurate to within .3°K and 
temperatures in the region of 160°K accurate to perhaps 1° or 2°K. A plot 
of these points for all three nozzle geometries tested igs shown in 
Pieure 31. 

Several characteristics of the condensation profiles and the local 
conditions at point of occurrence may be noted. 

First, the shape of the profiles contain no sharp minimum at the 
point of occurrence of condensation such as those noted in the condensation 
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of pure water vapor. .€. it 1s not appropriate to speak of condensa-— 
tion shock in this case, since the formation of moisture occurs along an 
extended region. 

Second it appears that pressure level (or temperature level, or some 
dependent property) may influence the shape of the profile, i.e. the rate 
of formation of moisture, to some extent since the shapes of the condensa- 
tion profiles appear to be more "Spread out” for high stagnation pressures 
than for low stagnation pressures. See Figure 27. These effects are not 
explainable in terms of gas dynamics influences, and the cause is not 
apparent. 

Third, a very definite influence of nozzle geometry upon the point of 
onset of condensation occurs, as seen in Figure 31. The rapid rate of 
nozzle divergence and expansion of Nozzle I results in higher supersatura- 
tion ratios and relatively greater delay of condensation. Progressively 


decreasing the rate of expansion through decrease of nozzle divergence 
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condensation (verified visually in Nozzle I). By noting the pressure ratio 
at which condensation occurs, it 1s possible to determine the local condi- 
tions at occurrence using the isentropic relations for the equation of 
state previously discussed. Based on the range of applicability of Plank's 
Equation of State and comparison with the tables of Din, it is estimated 
that temperatures in the vicinity of 190°K are accurate to within .3°K and 
temperatures in the region of 160°K accurate to perhaps 1° or 2°K. A plot 
of these points for all three nozzle geometries tested is shown in 
Figure 31. 

Several characteristics of the condensation profiles and the local 
conditions at point of occurrence may be noted. 

First, the shape of the profiles contain no sharp minimum at the 
point of occurrence of condensation such as those noted in the condensation 
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of pure water vapor. .e. it is not appropriate to speak of condensa- 
tion shock in this case, since the formation of moisture occurs along an 
extended region. 

Second it appears that pressure level (or temperature level, or some 
dependent property) may influence the shape of the profile, i.e. the rate 
of formation of moisture, to some extent since the shapes of the condensa- 
tion profiles appear to be more "spread out” for high stagnation pressures 
than for low stagnation pressures. See Figure 27. These effects are not 
explainable in terms of gas dynamics influences, and the cause is not 
apparent. 

Third, a very definite influence of nozzle geometry upon the point of 
onset of condensation occurs, as seen in Figure 31. The rapid rate of 
nozzle divergence and expansion of Nozzle I results in higher supersatura- 
tion ratios and relatively greater delay of condensation. Progressively 


decreasing the rate of expansion through decrease of nozzle divergence 
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condensation (verified visually in Nozzle I). By noting the pressure ratio 
at which condensation occurs, it 1S possible to determine the local condi- 
tions at occurrence using the isentropic relations for the equation of 
state previously discussed. Based on the range of applicability of Plank's 
Equation of State and comparison with the tables of Dingee it is estimated 
that temperatures in the vicinity of 190°K are accurate to within .3°K and 
temperatures in the region of 160°K accurate to perhaps 1° or 2°K. A plot 
of these points for all three nozzle geometries tested igs shown in 
Hieure 31). 

Several characteristics of the condensation profiles and the local 
conditions at point of occurrence may be noted. 

First, the shape of the profiles contain no sharp minimum at the 
point of occurrence of condensation such as those noted in the condensation 
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of pure water vapor. .e. it is not appropriate to speak of condensa- 
tion shock in this case, since the formation of moisture occurs along an 
extended region. 

Second it appears that pressure level (or temperature level, or some 
dependent property) may influence the shape of the profile, i.e. the rate 
of formation of moisture, to some extent since the shapes of the condensa- 
tion profiles appear to be more "spread out" for high stagnation pressures 
than for low stagnation pressures. See Figure 2{. These effects are not 
explainable in terms of gas dynamics influences, and the cause is not 
apparent. 

Third, a very definite influence of nozzle geometry upon the point of 
onset of condensation occurs, as seen in Figure 31. The rapid rate of 
nozzle divergence and expansion of Nozzle I results in higher supersatura- 
tion ratios and relatively greater delay of condensation. Progressively 


decreasing the rate of expansion through decrease of nozzle divergence 
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results in earlier condensation. Comparative geometry and non-condensing 
pressure profiles of the three nozzles may be seen in Figure 32. 

Fourth, the extremely low rate of divergence of Nozzle II in the 
0.80° case (subsonic flow results at about 0.50° divergence) results in 
notable axial fluctuation of the pressure profile data (Figure 29) and 
some scatter of the point of onset of occurrence of condensation. (Figure 
31) It should be remembered that differential changes in gas properties 
are proportional to (1/1-M*) and that at the pressure ratios present in 
these tests Mach Number remains very close to unity. (M = 1.25 at 
Be = .40). Also, as heat of condensation is released, Mach Number is 
driven even closer to unity. This is coupled with the fact that at low 
angles of divergence, the value of the negative pressure gradient is 
greatly decreased, downstream of the throat, from those values which 
existed for the other two test nozzles. This would tend to create more 
rapid boundary layer growth and greater departure from l1-dimensional flow 
in the narrow angle case. It might be expnected that larger departures 
from 1-dimensional flow, plus the magnifying effects of near unity Mach 
Numbers could contribute to larger pressure fluctuations. 

Cae Visual Observations 

Experimental observations indicated that the pressure hump accompany- 
ie cOndensavi eon Yoccurred dt~4. point COInci dent with Ehe appearance or 
the ice cloud. 

This at first seemed at variance with the expectation that pressure 
variation should occur prior to the appearance of visible condensate, 
due to the theoretically small sizes of forming nuclei. (r Pear ny. 
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investigation of condensing steam was also accompanied 
by visual observation of the water cloud coincident with the location of 


onset of condensation as determined by pressure measurements. He 
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reported that this observation remained valid despite changes which resulted 
in the axial shift of condensation in the nozzle. The same was true in this 
investigation. 

It is also mentioned that the data illustrated a very high degree 
of consistency and reproduceability, that illustrated in Figures 24 
through 28 involving three separate sets of measurements over a period of 
six months, each set being separated by nozzle dismantlement, cleaning, and 


pressure tap plate alterations. 
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IV. INTERPRETATION OF EXPERIMENTAL RESULTS 

As discussed in sections 1 and If, there are reasons for supposing that 
the process of nozzle condensation is one involving homogeneous nucleation 
within the super-saturated vapor. The purpose of this section is to discuss 
homogeneous nucleation theory and the subsequent drop growth process and 
apply these processes to nozzle expansions for the purpose of testing the 
ability of existing theory to interpret the experimental results of CO, 
condensation. Attention is given to the physical properties upon which 
nucleation and drop growth theory depend, and to the uncertainties which 
exist in the applicable values of these properties as well as to the theory 
itself. 

A. Nucleation Theory 

A review of classical nucleation theory, following the development of 
seein is presented in Appendix B. For more comprehensive review of 
theory see Reference Numbers 5,6,12,18,38,56,58,65 and 68. The discussion 
of vanes accomplished while working under P.P. Wegener at Yale University, 
is particularly thorough. 

LF General Form of Nucleation Rate Equation 

The rate of production of nuclei in a supersaturated vapor has been 

shown by various authors to be of the form* 


J = Aeexp) ~SO*/KT) (>) 


where J = number of nuclei produced per unit time per unit volume, 

AG* = free energy of formation of a critical cluster (a cluster in 
unstable equilibrium with its supersaturated vapor) 

And A is proportional to the incident flux of vapor molecules on a 


critical-sized cluster. 
- —==— x2 
A Cn ae ] unr (6) 


* See Appendtx B for a review of nucleation theory development 
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The value of C varies somewhat as derived in the classical treatments 
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Ot: veimerra: Becker and Deranene Frenkel, Bane Yang and others. 
Though variations of one to two orders of magnitude occur over the range of 
values presented, the effect on nucleation rate is small,since for only a 
small additional expansion along a nozzle, the value of J can easily vary by 
5 - 10 orders of magnitude. Since the change in local conditions is small 
for a variation in nucleation rate of this sort and, in fact, is scarcely 
noticeable for 1 to 2 order of magnitude variations in J, the differences 
in classical treatments of the theory are not significant. 

In the classical treatments of the theory, the nucleus is consi- 


dered as a stationery object, and the free energy of formation is merely: 


Je 
(GG = hnr2o - — 
A Nre=o eKT ln p (7) 


00 
where the first term is the work of formation of surface area for a 
liquid having surface tension o, and the second is the constant-tempera- 
ture free energy change of g molecules which change state from vapor at 
Eel to. saturated. liquid at. i. Ee is the flat film vapor pressure of the 
liquid at temperature T. 


- 
The -errvencall anise is that for which the function AG is a maximun, 


and results in 


Using, the procedure presented by Frenkel, for which the value of c 


becomes 


ce [paglbe yi? (9) 


the expression for nucleation rate is 
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As mentioned, Yang has shown that the results of the authors previously 
mentioned yield expressions for classical nucleation that are much alike 
(in some cases identical) and which yield quantitative predictions that are 
mou significantly different. 

More recent authors have pointed out the need for inclusion of addi- 
tional terms to the free energy of formation of a critical cluster, which 
account for the fact that the cluster may have translational, rotational 
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and vibrational energy modes. Kuhrt, Rodebush Lothe and Pound, and 


2st Shave treated this problem but reach widely different con- 


Courtney 
clusions as to the effect upon nucleation rate. The postulated correction 
factors for this "gassification" of the condensed clusters range from 104 
to TO. depending on the treatment and molecular configuration chosen. 

More recently, Feder, Russell, Lothe and Pounal® have jointly made a 
very comprehensive review of this problem, and offered some new thoughts on 
the subject. For water vapor, they conclude that the effects of gassifica- 
tion are equivalent to a correction factor of about ego but admit that 
this is not borne out by cloud chamber data. In general,such experimental 
data as does exist tends to more strongly verify classical eeeorn aoe 
than the improved gassified theory, despite strong arguments for the greater 
theoretical validity of the latter. 

fale Evaluation of Classical Nucleation Rate Equation 

Essentially the same physical properties appear both in the exponen- 
tial and pre-exponential terms of the nucleation rate equation. Of course, 
uncertainties in these properties will effect the calculated value of 


nucleation rate chiefly through the exponential. By substitution for the 


value of r* it is seen that the exponential is 
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exp! 3KT = exp!3 ‘GInP/P,’ (op) (jas) 
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In evaluating this term correctly it is clear that four physical 
properties of the vapor must be known: 1)the liquid drop surface tension 
2)local temperature 3)liquid drop density and 4)local supersaturation 
FAC1LO. 

Much attention has been given to the applicable values of surface 
energy for small droplets. However, here again there is little agreement 
as to what corrections should be made to flat film values to take account 
of the effect of droplet curvature, and the sign of the correction is itself 
a subject of controversy. 

Effects of curvature are considered by some to cause a surface tension 


decreasing with drop radius such as in the Kirkwood-Buff Equation. 


Do 


Gear) (a2) 
where 6 is a length, nearly constant, lying between .25 and 0.6 of the 
molecular radius. 
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In addition, Stever and Rathbun have calculated a correction to take 
account of the radius dependency of surface tension which also results ina 
decrease for small radii. Compared with uncorrected surface tension these 
corrections lead to lower energies of formation of nuclei and hence higher 
nucleation rates for the same supersaturated conditions. 

However, Oriani and Suid ietee predict an increase in surface energy 
of about 25% in the case of water, using an analysis based on the breaking 
of intermolecular bonds. A bonding model used by Benson and Shuttleworth, 
however, predicts a 15% reduction in surface energy of correspondingly small 
water ice crystals. Bogdonoff and acon also using considerations of bind- 
ing energies of sodium and of nitrogen obtained a correction that results 


in larger critical clusters and also larger energies of formation, hence 


lower nucleation rates from the uncorrected case. For lack of a more 
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definite indication of the effects of small droplet sizes, most investiga- 
tors use the flat film value of surface energy, and estimate the effects 
on the theoretical predictions of departure from it. 

An additional uncertainty in surface energy occurs with CO, due to 
the fact that, at pressures below its ice point, 75.1 psia, it apparently 
condenses in the form of ice crystals rather than liquid droplets. Bondi? 
has shown that for metal crystals, which involve valence bonding, increases 
in surface energy from the liquid value are closely proportional to the 
heat of fusion and, depending also on the crystal planes involved, amount 
to increases of about 25% for those metals which he discusses. He further 
presumes that for molecular solids one might expect increases of from 
approximately 1/8 to 1/3 these amounts. 

In the case of some gases there is doubt as to the proper value of 
droplet density to be applied. Though density does not vary too rapidly 
with temperature, there is doubt, in the case of CO,> as to whether solid 


or liquid density should be used. Although the CO, may condense as ice 


2 
particles, there is some argument for use of liquid density in order to be 
consistent with the spherical geometry and liquid surface energy applied to 


the liquid drop model, and in the absence of any knowledge of CO, crystal 


a 
surface energies. Also, studies of Dorsch and Hacker?© with very small 
water droplets indicated that the spontaneous freezing temperature 

decreases rapidly with decreasing droplet size. (Figure 33) This might tend 


to justify doubt on the reality of attempting to distinguish the difference 


between liquid and solid states at droplet sizes of critical clusters. [For 


CO,» typical critical radii ‘are of (.5°x 107° em to 107! em (about 30 - 50 
molecules) at supersaturation ratios of the order of 10, where noticeable 
a, 


quantities of moisture are condensed. Typical nucleation rates are 10 
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Figure 34 illustrates existing data on vapor pressures of CO, in the region 


a 
of experimental condensation, and the computational approximations used in 
the condensation calculations. Accuracy of the vapor pressure data and the 
approximation below temperatures of about 180°K is uncertain, but is probably 
within 5% for higher temperatures. Supersaturation ratio may be determined 
within these limits of accuracy for a given local pressure if the tempera- 
ture is also known. 

Surface tension is considered to be described for flat film liquid by 
the function o, = 5.00 = .193T(°C) dynes./cem, based on data of References 21, 
27, and 50, for which a linear dependence of gw on T down to 220°K is 


)++">( Reference 27) yields 


justified. An alternative function o = 75 (1a 
values of ow which vary by only 2% to 9%, in ae 200°K to 160°K region of 
interest, from those of the linear function above. 

Knowledge of local temperature, previously discussed, is important in 
evaluating the nucleation rate equation, not only because of its presence as a 
eubte ti the oxmonent, Dub because the other variables are also dependent in 
value upon temperature. Effects of uncertainty in temperature do not cancel 
appreciably in their influence upon theoretical nucleation rate. At 200°K, 

a 1% error in temperature will result in about a 2% error of the opposite 
sign in surface tension and hence about a 3% error in the ratio o/T for 


CO These effects were illustrated in the earlier discussion of the 


a” 
potential uncertainty of local temperature unless an equation of state is 


sea Co determine the 1sentrope for CU It is possible that some past 


O° 
applications of theory to experimental condensation data have erred through 
failure to consider these real gas effects. 

Departures from perfect gas behavior affect the evaluation of the free 


energy required for the formation of a critical cluster, and also the drop 


growth relations dealt with in the following section, since both nucleation 
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and drop growth developments are based on perfect gas assumptions. This 
problem is dealt with in Section IVD, following the discussion of drop 

growth theory, and the application of nucleation and drop growth to condensa- 
vaon in nozzles. 

Be Drop Growth Theory 

Subsequent to the production of nuclei, additional formation of 
moisture will occur through the growth of the clusters due to vapor impinge- 
ment. Nuclei may by this process grow to radii of up to about an order of 
magnitude greater than their original size, depending upon two important 
drop growth parameters introduced below. In the theoretical treatment, the 
distinction between nucleation and commencement of drop growth treatment is 
somewhat artificial, as it is necessary to establish some size of cluster 
at which the nucleated droplet exists and may commence growing. Becker 
and Doring have shown that at a cluster radius 1.3 times critical, the 
probability of clusters continuing to grow reaches better than 99%, and 
hence it is considered that this is a reasonable starting point for the 
application of drop growth equations. (Note: Computations showed that 
results were relatively insensitive to this starting radius.) 

The growth of a droplet depends predominantly on two factors: 1)The 
rate of condensing mass flux 2)The rate of droplet energy transfer to. the 
vapor environment. 

It would appear that the first factor is itself a description of the 
drop growth process, but in actuality, it is greatly dependent upon not 
only local vapor conditions, but also the difference between droplet and 
vapor conditions. Both of these are strongly influenced by the second 


factor. 
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Two coefficients are defined to describe this behavior: 


& = Condensation Coefficient. The fraction of impinging vapor 
molecules which are absorbed at the drop surface. 
a = Thermal Accommodation Coefficient. The fractional temperature 


change which takes place in the reflected portion of incident 





molecules. 
eee s 

ee (13) 
T- sh 


Temperature variation across the droplet is assumed to be negligible, 
and the droplet is assumed to be spherical. 

Based on mass and energy balance considerations presented and developed 
in Appendix C the following drop growth equations are obtained. 


Conservation of Mass 





P 
DOr is 2,___EP Reece (14) 
Pas > TE aya /2 (ant yt/2) 
Conservation of Energy 
aT 
yey sea) 2dr = 7 
4/3nr p,C ee ae Lar at ep B( ORT) EB (QRT. ) 
- (1 -&)82R(T + a(T, - T)) (15) 
These equations may be reduced, respectively, to 
chou 
oa (16) 
Z BD Tp Tp Boy _¥ 
OS PAN TM rs ape) AOL ES Ae ee SS Ne peg a) eee 
a X= rhe 
where A = z RT 
Jo gs a oe 
= (ompp)l/2 De (27RT,) He 


And Py is the equilibrium vapor pressure for a spherical liquid drop of 


temperature T, in an environment of temperature Th: 
OO 
7 (————) 
oe a) exp 0, RT Lr (18) 


In the further application of these equations to streamwise growth of 


drops in a nozzle, additional assumptions are made. First, for an 
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acceptably simple computation scheme it is convenient to treat droplets of 
uniform size, where the radius i8 an average based either on surface area 
or volume. Ignoring the effects of coagulation, there is evidence that the 
actual distribution of drop sizes is narrow, and that the average droplet 
treatment does not introduce appreciable error. 

One possible error that may occur as a result of the above assumption 
becomes likely when the supersaturation ratio P/P,, is reduced so rapidly 
as a result of condensation, that droplets, which have recently nucleated 
and begun to grow, will suddenly be of size less than the critical cluster 
size. In this case, considerable re~-evaporation would occur among the 
smaller droplets of the actual size distribution, and the treatment of an 
"average" droplet does not account for this. 

Relative motion between droplets and the gas stream is ignored, and 
though this is generally accepted, it may be possible that the same 
translational modes of energy required by recent authors in nucleation 
theory might also affect to some degree such things as incident mass and 
energy fluxes upon newly formed nuclei. 

Values of q and &€ have been presented by numerous investigators for 


several vapors, excluding CO Although earlier presented values of & were 


elt 
fairly low, there has been more recent substantiation of values close to 
1.00 for most materials, including liquid metals, the previous errors 
having been caused by neglect of surface temperature depressions and the 
effects upon evaporation flux. However, there is also indication that 
molecules of higher dipole moment tend to have lower values of &, as with 
eater ; 2 9 
liquid water for which Alty and MacKay, and Prueger present values of 


less than €= .04,. The influence of crystallization upon & may also be 


large, as Tschudin presents a value close to unity for water ice. Effects 
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of surface contamination and surface diffusion barriers (to recently 
evaporated particles due to insufficient vacuum and mean free path lengths,) 
are also thought to play a role in experimentation which tends to lower 
apparent values of € No experimental values of & have been measured for 
CO,» which has a dipole moment of about 1.5 compared with 80.0 for that of 
water. 

The thermal accommodation coefficient for vapors and their liquids 
is generally accepted to be around unity. However, values may be much 
lower between vapors and solids, as Poems indicates for gases on clean 


metal surfaces. For CO the value of a for the vapor and its solid 


0°? 
crystal is unknown, and there is uncertainty as to which of the two situa- 
tions above is the most applicable description. 

Departures from perfect gas behavior also influence drop growth. 
Again, correct values of local temperature and supersaturation ratio will 
be influenced, and the perfect gas assumptions used in treating mass and 
energy balances leading to equations (16) and (17) may require modification. 

In this study, considerable attention has been given to the behavior 
and influence of the above mentioned factors. 

on Application of Nucleation and Drop Growth Equations 

The predicted effects of moisture formation during a nozzle expansion 
are determined by applying the time rate release of heat to the gas dynamics 
of the related expansion. By considering the process with respect to 
variable nozzle position, rather than time, a result is obtained which can 
be directly compared with experiment. 

Interferometric measurements of density variations in the nozzle 
verified that the flow between boundary layers is in fact one-dimensional 


and allows the use of the one-dimensional gas dynamics relations; 
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es eT (19) 
Momentum 

-Adp = mdu (20) 
Energy 

C tat = udu = vh. du (2m) 

p fg 


where m is the mass rate of flow through the nozzle, H the mass fraction 
of condensed moisture, u is the streamwise velocity of the vapor, and 0, 
P, and T are the vapor density, pressure and temperature, respectively. 
A is the effective area of flow between boundary layer thickness, and is 
determined experimentally for non-condensing flow. 

The rate of moisture formation is determined as follows and assumes 
the droplet velocities to be those of the vapor. A nucleus of radius re 
formed at some point Xo will grow to a larger size at some future time and 


adistance x. 


and the surface area will be 


(é 
= Ante + x = — ax]? (23) 
dr 


For a rate of growth at 


predicted by drop growth, the increase in 


liquid mass per droplet will be 


dr axj2 22 i 
: dx | ae dx (2h) 


ome oe 


x 
adr = 9, dol. + l 
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And the fractional moisture change that occurs at x for the 


dn = J(x_)A(x_)[——~ ]dx_ particles formed in the volume A(x _)dx_ is 
O O Grau O SO 
° aXo x j dr ba es 
= dnoPLadr _ J(Xo)A(Xo) Uo axoep,tlrot lx, u_dt dt_u 
a a ae rae = 
O P A dx,)/ He uo? (p Agdx,)/(2 u,) 
(25) 


Inclusion of the drop growth on all previously formed nuclei necessi- 


tates integrating the right hand term over all Xo° Noting “that 
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nm = © AgY/t - uy and considering also the moisture produced through 


additional nucleation at x, there results 


Cie” 3: A( Xo) dX9)0 L4n | pal eche 2 dr 
ax I(x) ( ) m Ir, y xu dat dx | at cee 
Sai (8) 

+ HU AL ox) A(x) 

3 lt 
_ WvPy, 1 dar > ede 
oh UY O I dt oh) INGA u dt —%o 
+ 1/3 r° J(x)A(x)] (26) 


For nuclei generally only slightly larger than critical size, the 


theoretical drop growth rate is essentially independent of ete Hence 


dar 


at is assumed to have a value equal to that for large drops. Also, since 


it 1s practically impossible to keep account of the radii of each group 
of growing nuclei which formed at each section dx, it is assumed that the 
droplets can be represented by the "surface-area averaged" radius of 
droplets. As previously mentioned, volume averaging has little influence 
on the results. In addition, Lundgren has given experimental evidence 
that the actual size distribution of particles condensed in a nozzle 
expansion is quite narrow. 

The foregoing equations are reduced to the following form for use 
with a Runge-Kutta=-Merson integration procedure along the expansion path. 


(The detailed listing of computation procedure may be seen in Appendix H). 
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where 
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Initial values of the stagnation conditions, physical properties of the 
vapor and its liquid, values of a and €, and calculation starting point are 
given for each computation. Effects of contamination may also be considered 
by giving non-zero values to Y(1), Y(2) and Y(3), based on an initial concen- 
tration and size of contaminant. Effects of possible droplet energy modes 
upon free energy of cluster formation ("gassified" behavior of droplets - 
page 23) may be considered by multiplying nucleation rate by constant arbit- 
rary factors, equivalent to the values proposed by various authors. 

The various sub-routines of the listed computation scheme accomplish 
the determination of such things as the change in nozzle area, vapor and 
liquid properties, supersaturation ratios and solution to the drop growth 
equations when step-wise changes are made along the stream of the expansion. 
An error estimation is made for each computation and the size of the step, 
dx, is adjusted to reduce the error bounds to within the desired limits. 

Ds Real Gas Effects 

The preceeding developments of nucleation and drop growth theory and 
application to nozzle condensation computations assume perfect gas behavior 
during the expansion. However, as earlier discussion illustrated, experim- 
ental measurements were being made for which values of =~ which were consid- 
erably less than unity. Consequently, the computational scheme was modified 
so that the effects of departure from perfect gas behavior could also be 
considered. 

A feature of an isentropic expansion which facilitates fairly simple 


corrections to the gas dynamics, drop growth and nucleation equations is the 
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fact that for this particular thermodynamic process, the compressibility 

factor remains fairly constant, changes of less than 2% in the value of z 
generally occurring for expansion pressure ratios of 10 or more. This was 
found to be true for several other gases, including ammonia, and for CO, 
is valid throughout the superheated and supersaturated regions for which 
compressibility factors are in excess of about .88. Using this fact, it 
may be shown (Appendix D) that the following relations are approximately 


true for an isentropic expansion. 


Pv 


Rr zZ =constant C210) 
_ Bae 

Y= Tete) * ao (28) 
Ri = <2 = gp (29) 


where y' and R' are the effective ratio of specific heats and gas 
constant, respectively, that are applicable. It should be noted here, that 


the nucleation rate is effected in the following manner; 
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due to the change in free energy required to form a critical cluster. 

To apply these effects to the perfect gas relations of the condensation 
computations, two additional steps were necessary. First, modifications which 
allowed for use of a variable ratio of specific heats (Y = Y(T)) were made. 
Second, each set of stagnation conditions which was to be placed in the 
condensation scheme first had its "real" isentrope, z, and effective isentro- 
pic exponent Y'(T) determined, and from these the Y(T) which corresponded 
to z and Y'. Computation input was Y(T) and correction (compressibility 
factor) z. This resulted in an isentrope, providing the correct values of 
P vs. T being employed in each condensation calculation, plus the use of 
effective values of y and R in the drop growth and nucleation calculations. 


Additional details of this procedure are presented in Appendix D. 
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ine Results of Application of Theory to Condensation of CO» 
i ie Comparison of Theory and Experimental Data 


Data shown in Table 14 and Fig. 31 illustrate the local conditions at 
the occurence of condensation for each of the experimental tests made. 
Local temperature was determined through computation of the isentrope for 
each test case, using the procedure discussed on page 13 based on Plank's 
Equation of State for CO». Reference to Figure 35 illustrates that .1% 
presence of moisture is a convenient and reasonable criteria for theoretical 
onset of condensation. Comparing the data and the predictions of theory 
involved two main comparisons; 1) The location of the onset of condensation, 
defined as the theoretical occurrence of .1% moisture, and 2) The shape of 
the pressure profiles. For lack of better knowledge, %, §& and = were 
originally assumed to be unity. Specific comparisons were made of the pred- 
ictions of 1) Classical nucleation theory 2) Droplet "gassification" 
corrections to the theory 3) "Real" gas vs. perfect gas treatment of the 
classical theory (to determine whether or not a simplified approach, contain- 
ing known errors, might yet yield a better prediction of observed behavior. ) 
4) Use of liquid vs. solid droplet density and 5) the predictions of the 
theory in nozzles of different geometry. Figures 36 through 39 illustrate 
these comparisons and support the following conclusions; 

1) All predictions based on values of 5/0,, 4 and € of unity 
predict condensation at considerably lower supersaturation ratios than those 
observed experimentally. 

2) Real Gas treatment of the classical theory yields a closer 
prediction than perfect gas treatment, and the "gassification" corrections 
are in the wrong direction. (Fig. 36) 

3) Use of liquid rather than solid droplet density yields a 


prediction nearer to experimental results, lacking any better knowledge of 
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what corrections might be applied to surface energies of small CO. crystals. 
(Fig. 37) 

4) The unmodified theory does extremely poorly in predicting any 
effects of nozzle geometry. (Fig. 38) 

5) The shape of condensation profiles predicted by values of @ 
and © of unity are in extreme disagreement with those observed experimentally. 
(Figs. 45 and 46) 

2. Influence of Nucleation and Drop Growth Theory on 

TWaeoreticals Prediecions., 

kt is concluded that homogeneous nucleation and drop growth theory fail 
in their ability to predict observed condensation behavior of CO,» unless 
some modification in values of the nucleation and drop growth parameters oa, 
Pros a, and € can be independently established. Toward better understanding 
of how these parameters influence predicted condensation phenomena, and with 
the idea in mind of finding corrections for the prediction of CO, condensa- 
tion, a series of calculations based on arbitrary variations of a, Pr» o and 
— were made. 

Figures 39 and 40 show typical pressure profiles for condensation of 
CO,» for various condensation and thermal accommodation coefficients. Shown 
in the lower portion of the figures is the nucleation rate, indicating how 
nucleation rate remains high for longer duration when drop growth is retarded. 
This results in condensation clouds of greater particle density and smaller 
particle size than in the case of ready drop growth. In addition, smaller 
total deflection in the static pressures (from non-condensing), indicating 
lesser amounts of moisture, exist for the downstream portions of profiles 
having lower drop growth parameters. 


An important observation in these figures is that, beyond a certain 


reduction in the drop growth parameters, appreciable proportions of moisture 


“a 
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may be present primarily from nucleation. This is manifested in the form 
of a so-called "ramp" deflection in the pressure profile to which additional 
drop growth effects are super-imposed. Thus, for example in Figure 39, for 
an assumed a of .10, variations in & from 1.0 to .3 and then to .1 have 
progressively less effect on the point of departure of the condensation curve, 
and for lower values of &, the nucleation "ramp" becomes increasingly appar- 
ent. The same behavior occurs for decreasing values of a, (shown for an 
assumed value of &€ = .04) in Figure 40. It would be possible to obtain a 
deflection curve involving only nucleation, if drop growth were set at zero. 
Figure 41 illustrates the influence of increasing or decreasing the 
nucleation rate by changing the free energy of formation of nuclei. Indicated 
are variations in the surface energy from flat film liquid value, though 


Similar changes may be introduced by varying the droplet density 0 The 


L° 
important parameter in the nucleation rate exponent is (o7/oee)) since, with 
knowledge of temperature and supersaturation ratio, the uncertainties of 
surface energy corrections, small cluster densities, and non~spherical geom- 
etry may all be contained in this term. 

Two interesting features may be noted in Figure 41. First is the effects 
of higher downstream velocities, which tend to "wash out" the curves further 
downstream, and second, the important observation that, for a and € of unity, 
condensation results in a return to near-equilibrium conditions, seen by the 
fact that condensation commencing at any point in the nozzle results in the 
same pressure deflection immediately further downstream. 

a) Onset of Condensation 

As the previous figures suggest, influence of the drop growth coeffic- 

ients upon the point of onset of condensation is relatively small compared 


to that caused by changes in surface tension. However, uncertainties in a 
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and & still have appreciable effect upon values of surface energy necessary 
to bring the predictions of theory into agreement with experiment. Figures 
nO, 43, and 44 show the corrections to flat film surface tension necessary 
40 provide agreement between theory and experiment (based on fitting the 
data of nozzle IX2.20°)). i the uncertainty in the values of o and & for 
CO, crystals, is given limits of from perhaps .01 to 1.0 for each, then the 
respective surface energy corrections necessary may range from about 20% 
to 38% increases in the flat film surface tension. For the absolute lower 
bound of uncertainty on a and € of zero, the necessary correction in surface 
energy is reduced to +17%. 

kt is noted in these figures that the theoretical effects of nozzle 
geometry are increased to some extent by use of small, but non-zero values 
of the drop growth parameters, but that the theory in all cases predicts less 
effect than that recorded experimentally. 

b) Profile Shapes 

Figure 45 illustrates the pressure profiles for a single exparsion 
corresponding to the same values of c@ and €, and corresponding o determined 
above, and indicates an additional distinct failing of the theory in its 
description of the early stages of condensation. kis apparently impossible 
to reconcile the experimental downstream pressure deflections, which indicate 
a considerable degree of supersaturation, or departure from equilibrium, with 
the apparent rates of formation of moisture immediately after onset. When 
drop growth parameters are adjusted to bring theoretical profiles into 
agreement regarding the amount of downstream pressure deflection, then the 
upstream theoretical rate of formation of moisture lags considerably behind 
the experimental results. 

This problem may result from some of the approximations discussed in 


the numerical calculations, or may in fact indicate some physical "arrest” 
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of the condensation process, perhaps due to decreasing mass or thermal 
accommodation characteristics of crystals which grow rapidly colder during 
their passage downstream. 
c) Downstream Moisture Formation 
For purposes of calculating the amounts of carbon dioxide that will 
form, or the downstream pressure variations which will occur during a nozzle 


expansion of CO,, the following parameters were determined, having arbitrar- 


y iad 
ily set thermal accommodation coefficient, a, at unity. There results, as 
a best description of the three nozzles, the profiles shown in Figures 46, 


47, 48 and 49 for a = 1.0, & = .00008, p, = 1.18 gm/em?, p/p, = 1.18. These 


L 
drop growth parameters have no physical significance, and are presented only 
as an aid in determining approximate amounts of CO, moisture formation sub- 
sequent to the occurrence of condensation in a nozzle expansion. 

The best compromise to profile shape soon after onset has been a recipe 
which yields theoretical profiles lagging behind and below the experimental 
data early in the condensation, and which then pass slightly to the other 
side of the data by about the same amount for the remainder of the profile. 

The condensation data of Nozzle II (0.80) indicates the presence of so 
some cyclic axial variation in the pressure ratios of all profiles. This 
suggests that at this narrow angle, some departure from 2-dimensional flow 
may be interfering with the condensation process. Friction, or compression 
shocks may be creating these fluctuations, and would also result in decreased 
supersaturation ratios in this region, offering a possible explanation as to 


why considerably less moisture is forming than the recipe predicts for that 


nozzle. 
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3. Influence of Contamination 
Computations attempting to determine the influence of contamination of 
vapor in nozzles were carried out at intervals during this investigation. 
Original calculations were based on perfect gas assumptions and values of a, 


E 


2? 


sand = of unity, in an attempt to establish whether or not condensation 
could ates be due to reasonable amounts of contamination present in the 
test vapor. Such contaminants might be present in the form of dust, charged 
particles, or the result of an earlier condensation of some other vapor. 
Contamination of this sort would be expected to cause earlier condensation 
than would occur if the test vapor were perfectly homogeneous. 

none presents data for typical size density distributions of contam- 
ination in atmosphere. Upper limits of particle concentration in dirty 


“3 for particles of 


atmosphere of a large city reach values of 5 x 10° em 
size less than 10° em (.1 micron), about 10° cn? for sizes from .1 to 1.0 
microns and 1.0 cm? for particles exceeding one micron in size. 

Typical profiles of supersaturation obtained for increasing concen- 
trations of contamination with microscopic particles are shown in Fig. 50. 
Beyond a limiting concentration of particles, (depending on size) the super- 
saturation ratio obtains lesser and lesser values, until the flow reaches 
a state where it is close to equilibrium throughout the expansion, due to 
the ready formation of moisture on contaminant particles. 

Subsequent to the establishment of the "recipe’ described in the 
previous section, additional computations for both microscopic and macro- 
scopic contamination were made. 

The results are shown graphically in Figures 5l and 52, and are briefly 
summarized: 


3 


1) Micron-sized particles in excess of 108/ em? are required to 


influence the theoretical condensation of CO.,. This is at least 5 orders of 
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magnitude larger than upper atmospheric concentrations of particles of this 
$size. 

2) Microscopic particles, such as might occur from the preconden- 
sation of another vapor, of size the order of Omer are required in excess 


of 1012 /em?3 


to influence theoretical condensation of CO,. 

3) No appreciable difference between empirically fitted real and 
perfect gas theories results. 

4) The difference between amounts of contamination necessary to 
influence condensation in the rapid expansion of Nozzle I and the slow 
expansion of Nozzle II (0.80°) is approximately one order of magnitude, 
independent of particle size. 

5) The empirically fitted recipe predicts necessary concentrations 


of particles roughly two orders of magnitude larger than that predicted using 


values of &, &, and 9/0, of unity. 
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two curves tend to move together, since the first curve to drop downward 


When the curves of Fig 52 are replotted in terms of {log , the 
represents larger volume contamination. This results in fairly localized 
values of supersaturation ratio at occurrence of condensation as a function 
of volume contamination, regardless of particle size. 

Applying the volume amounts of contamination of CO, in nitrogen presen- 


| )y 
ted by Willmarth and Neeanet ene. and Arthur 3 


» a computation of the predicted 


effects of theoretical influence of contamination on nitrogen was made, using 
values of a, &, and a/o of unity and perfect gas assumptions for nitrogen. 
The results are shown in Figure 53, for various size particles typical of 
cluster sizes for CO. > indicating reasonably good agreement with their data. 
This is taken to indicate some substantiation of the above calculations for 
contamination in CO,, perhaps within one or two orders of magnitude of the 


calculated contamination levels needed to influence homogeneous condensation. 





Me SUMMARY OF FINDINGS AND CONCLUSIONS 

A. Experimental Findings 
1. Consistent data regarding the supersaturation of CO, at onset of conden- 
sation plus the shape of pressure profiles has been obtained for three 
different nozzles. 
2. Interferogram photographs and an approximate theoretical estimation 
show agreement that the boundary layer in the nozzle throat was of the 
order .006" and substantiates the treatment of flow in the mid-stream as 
l-dimensional isentropic flow. 


3. Departures from perfect gas behavior by CO, in the P-T range tested are 


2 
Significant, and must be accounted for in the calculation of local conditions 
along an isentropic expansion. Density and static pressure measurements 
plus total pressure measurements for non-condensing flow indicate that the 
Plank Equation of State gives a good description of low temperature CO, 
behavior. 
4, Rates of CO,, expansion, as determined by nozzle geometry have a large 
influence upon the degree of supersaturation obtained prior to occurrence 
of condensation. 

Be Application of Nucleation and Condensation Theory 
1. Classical nucleation and drop growth theory, as applied herein, will 
reasonably predict the conditions at which onset of CO, condensation occurs, 
if flat film liquid surface energies are corrected by amounts within 
suggested ranges of uncertainty for small droplets or crystals. In these 


caiculations a droplet density equal to that of CO, liquid at the triple 


2 
Point. 16 em/em™ , has been assumed. 
ec. The appropriate correction to surface energv is dependent upon the 


values of mass and thermal accommodation coefficients, € and a, and has a 


correction uncertainty of about 21% of flat film values for maximum possible 


uncertainty in & and a. G Se boi(cror s—(@ = 10L0. = “= ..sc tor 
oo oo 
— =a: 1.0.) Consequently, the applicable value of the parameter 
ig 70, >) to be applied to the nucleation rate equation corresponds to 
values of (07 /oF) (0, */03) which lie between the limits (1.17)? and 


(1.38) 9where o, = 1.18 gm./cm3 and o| is the linear temperature depen- 
Oo 


5.00 = .1931(°C) dyne/cm, 


dent function o. 
3. Pressure profiles obtained experimentally indicate that considerably 
less moisture is forming than that necessary to return the mixture of 

vapor and condensate to near equilibrium values, or zero supersaturation 
downstream of the onset of condensation. Values of € and a necessary to 
produce agreement between the theory and experiment in this respect require 
that one or both of these parameters be considerably less than what are 
thought to be minimal values (about .01, for each.) 

4, Even at their best empirical fit, classical nucleation and drop growth 
theory as applied herein do not seem to provide any means of describing 

the exact shape of the pressure profiles obtained experimentally or the 
extent of effects of geometry noted experimentally, These failings, in 
addition to that of B3, above, are thought to indicate inadequate or incorrect 
assumptions in the treatment of the drop growth calculations. 

5. Application of the so-called "gassification" corrections to classical 
nucleation theory results in corrections in the wrong direction, unless 
even larger corrections in the surface energy of formation of CO, ice 
crystals (than indicated above) can be established. Use of solid, rather 
than liquid, droplet density in the parameter (eos) further increases 

the necessary corrections to surface energies required to provide agreement 
with experiment. 


6. The influence of dirt particles or other foreign particles on the 


a 


ie tom 


nT i ‘ 





condensation of CO, » and probably most other gases,may be completely 

disregarded in supersonic nozzles, unless the particles result from the 
precondensation of another gas in the stream. Amounts of micron-sized 
particles necessary to influence CO, condensation are about five orders 


of magnitude in excess of upper atmospheric limits of concentration, and 


particles of the order of 20 a° are necessary in excess of about sae 
Particles/cn.” to influence homogeneous nucleation and condensation. 
Testing the condensation of nitrogen with pre-condensed particles of CO, 
in quantities mentioned by Willmarth and Nagamatsu, and Arthur, yields 
theoretical predictions which are in reasonable agreement with their 
experimental observations. 

7. Additional knowledge of thermal and mass accommodation coefficients 
are necessary to allow better testing of homogeneous nucleation and 
condensation theory. In addition, comprehensive testing of the theory 
will await some better establishment of surface energy corrections for 
small particles. Until this information becomes available and an estab- 
lished treatment of nucleation and drop growth can be applied with con- 


fidence to many vapors, "recipe" fitting may provide the best means for 


predicting nozzle condensation behavior. 
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VI. SUGGESTIONS FOR FURTHER STUDY 

One of the greatest difficulties in application of experimental 
results to a meaningful test of theoretical predictions results 
from the uncertainties in value of many of the important variables 
and parameters contained in nucleation and drop growth theory. Con- 
sequently, any methods or devices of experimentation which would 
decrease any of the uncertainties, or the effects of uncertainty, 


of such quantities as 0, 9p., a, &, local temperature and the gas 


L 
dynamics in general would be of great assistance. One means of 
proceeding in this direction is to contain the test gas within large 
quantities of another gas which is non-condensing in the region of 
interest of the test gas, and whose gas properties are well-known, 
preferably being as near those of a perfect gas as possible. An 
immediate benefit of this is to provide a means of accurately determin- 
ing local properties at all points of the isentrope with ease, since, 
the temperature of the test gas should be that of the ideal "carrier" 
gas, which dominates the gas dynamics. Second, the influence of 
thermal accommodation is greatly reduced due to the fact that the 

huge preponderance of incident molecules upon the surface of a formed 
droplet are those of "carrier" gas, and hence are reflected. The drop 
is consequently kept at or nearly at environmental temperature 
regardless of the value of a. Hence, there are advantages to be 

gained by testing gases for which these uncertainties exist (the 

huge preponderance of all gases) in a carrier stream of such gases 

as nitrogen, helium, or dry air. Several complications arise from 

this method, including the need for very precise measurement of mass 
fraction of test gas in the mixture, and more sensitive instrumentation 


to detect the lessened effects of condensation upon the gas properties. 


a 
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Highly steady flow systems would aid in obtaining the accuracy 
necessary. 

With regard to some of the uncertainties pertaining to CO,» ae 
would be of considerable interest to determine whether or not the 
experimental behavior of condensation is influenced by establish- 
ing liquid rather than ice formation of condensate particles. This 
could presumably be accomplished by raising pressure levels to a point 
where condensation occurs above the ice point, 75.1 psia. This 
would necessitate stagnation pressure in’the range of 175 - 250 psia, 
mass rates of flow and a source of CO, which “are oroporti onal ly 
larger, heavier construction of some important fixtures in the system, 
(such as the stagnation tank) and a means of pressure profile procure- 
ment other than manometer board. 

It is felt that interferometric means of study offer enough 
advantages, particularly with more toxic or reactive gases, to warrant 
further application to study of condensation in nozzles. The inability 
to distinguish interference lines in the upstream region of the nozzle 
encountered at high pressure levels in this investigation could 
possibly be solved by facilitating higher camera shutter speeds 
(1/150 sec., the maximum possible, was used in these studies) or by 
establishing wall temperatures which are close to stream stagnation 
temperature. Two-dimensional nozzles, where interferometric study 
could be successfully employed, would decrease the problems of 
instrumentation and would be readily adaptable to studies of light 
scattering and polarization, yielding experimental information 
regarding droplet sizes. Interferometric photos involving conden- 


sation, although not applied in this report, were obtained, and showed 


interference bands very clearly within the region of condensation, 





? 





NT. 


this@ugeion bein~ marked plainly by the decreased light intensity 
Heoculiing trom the presence of the ice clond. 

Light scattering or other techniques which yield infor- 
mation regarding the average size of the ice particles at their 
initial appearance and after subsequent downstream growth would be 
very useful. Since, as previously mentioned, there is some indication 
that nucleation is the major contributor to moisture in its early 


stages of formation for CO,, this information might be directly 


5? 
melated. to the critical Cluster -size of CO, « 
Additional investigation of crystal thermal and mass accommoda- 
tion coefficients would be of help in further establishing the 
apparent indication from this investigation that either one or both 
have very low values for crystalline CO,, at low temperatures. ‘The 
carrier pas technique previously mentioned could isolate the influence 


of a and might facilitate a better means of studying € and its role 


in condensation. 





™ 
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FIGURE 7, ASSEMBLY VIEW AND GEOMETRY OF NOZZLE [. 
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FIGURE 8, DETAIL OF NOZZLE ILI 
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(C) FLOW, SHOWING SHIFT IN INTERFERENCE BANDS 


FIGURE 11, INTERFEROMETRIC PHOTOGRAPHS FOR TEST 1 OF DENSITY PROFILES. 
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(B) SCALING PLATE (1.90"') 
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(C) FLOW, SHOWING SHIFT IN INTERFERENCE BANDS 


FIGURE 12, INTERFEROMETRIC PHOTOGRAPHS FOR TEST 2 OF DENSITY PROF LILES. 
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FIGURE 17, EXPERIMENTAL NORMAL SHOCK CHARACTERISTICS OF CO5. 
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FIGURE 18 , PREDICTED NORMAL SHOCK CHARACTERISTICS OF CO. - PLANK EQUATION 


OF STATE. 


NOTE: POINTS CORRESPOND TO NUMBERED EXPERIMENTAL Py? To? AND P, 
EXCEPTING POINT #13 








10 


sat ATM, 


PRESSURE 


A2T 


FIGURE 19, COMPARISON OF PREDICTED CO, ISENTROPES 
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FIGURE 20, COMPARISON OF PREDICTED COy ISENTROPES 
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FIGURE 21 , CONSTANT NUCLEATION RATE PROFILES OF CLASSICAL 
NUCLEATION THEORY FOR CO,. < = 1.00. 
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FIGURE 22, EXPERIMENTAL PRESSURE RATIO PROFILES. NOZZLE I, 
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FIGURE 23, EXPERIMENTAL PRESSURE RATIO PROFILES. NOZZLE I, 
PRESSURE PLATES 1 & 2 


o #86TEST # 2 
° TEST # 4 
A TEST # 6 


—— NON-CONDENSING(TEST # 8 AND TOTAL 
PRESSURE MEASUREMENTS) 


x > 
\ 
ee 
‘ 
. o 
\ 

ae 

x A 

\ 

\ 

‘ R 
\ 
\ 
‘ 
\ 
N A 
\ 
~ 
nee a 
™~ 
eae A 


STATION NUMBER 


10 ll rz 13 14 15 


136 


16 





A33 


°94S933 aanssard [e,0} JUenbesqns BOI pautwsazap SonTea paancddwTt sayeotpuyr, 


e°d Ge 
2°6 OZS0° eee & he 
(oye G6¢90° are Ce 
9°6 oie G29C° €°Q ae 
-- c°6 C°3 3020° Le Te 
-- ake, 9°3 Ond0° T°6 O2 
Be 6°9 €°6 6920° 1°6 61 
-- O29: ¢°6 Cole 1° 6 900° 6°6 QT 
~090° == 9¢66° 6°6 £990" C*OT 930° 5 OL ei 
2290° oo Gots 6°1T ZOT° Beara 90T° om oa 91 
<060° == Gat" + ° ST HT’ esol ¢ GQST° Orel GT 
ste <= 6QT° 6°6T 9QT° lace Q9RgT° ter HT 
eTQT° we 6L2° 4°62 G2c° 6°c¢e bees C° He eT 
secc° -- Qs ° S°On O6¢ ° o° Oh T6¢° C°Qh ral 
Ca is L6n° C°2S TES: S°09 O2S° O°S9 Tt 
898s° $° 89 AOT 
pele 0°O° Gio. T°SS 449° 4°92 ess aS?) OT 
GQ?’ as i A6 
S9d° 4°29 92° eB", 992° 9°06 o22° G° 46 6 
49Q° L°rOT AQ 
600° T o- 12 O00" T 9°SOT 000°T eet al 000° T Orca 9 
0°6¢T AL 
Sa/a (3H )a 84/2 (FH )d Bais (FH. )d 84/2 Colao “ON UOTRIYS 
dgn'gnt Fd ecient ee “Sho = 
‘ age? dS nf AoSh- = Lb qe 64 dL 
. 5 (PAzeTNOTED) | a 
etsd L°G2= d ets? 2°26 = 4 ZuuC°2e = 4 etsd 2°2Tt 2d 
Ty son ihe *g °ON uny “+ “ON UNY To -on ung 


€¢2°ON aInStTy azoxy RYE auNnssacg 


9 Tava 





A3); 


ae 
Vv — 
—,, 
Mi , at 
A Vv a 
x —< 
9 aes 
x Vv Sa 
8 ™ 
vw 
x ae 
] ™~ 
x a 
° x ba 
@ x 
x oN 
8 : \ 
x 
oN 
Cs) 
a \ 
SON 
cad 
° 
oO 
Gv SSIS. = 


91 ST 71 tT cl Il Ol 6 
YadWNN NOLLVLS 


7 +¢ LSAL Vv 


€ # LSSL a 


c # LSaL ° 


I # LSA. x 


JO WIONVY ‘1 SLVId FunSswd ‘LI FIZZON 


O79 Ot 


9812 = FONTOUAAIA 
“STIIAOUd OLLVA FANSSHNd TWINDIUTEKT “77 TUNOTA 


SHHONI X 


Wah 





A35 








LOU? O° Ce ego G*te Qte° e°ek Se” C°0¢" 6d ce 
Old" @-ce  60c- T°9¢ nQ2° h°to =192° O°HS 492° 7°2¢ ie t2 
¢ec° 9°¢2 S2¢22° Roce hid” 9°82 Oe’ G°9¢ s_  GZ2@* GE” O° He Ge Oe 
62¢° ¢°S2 ge’ Ce Oee CGc. €°0f ¢0¢° G°6E 662° O°6£ #£062° 5 °9¢ Ce 6T 
9Se° Co2d* -OS¢c° ¢:2o- 0Ge* O°Ss TEs? T°Sh 390 ee” C°eh 4TS" 6°6¢ Te che 
yer C766 WOR. O°S¢ coc’ Gk GE? G°9h TSE" L°Sh Onf' T°Sh 6T ck 
QTS° 9°¢e = QT” L°6¢ ¢oE° ¢°6E 236%" 6°TS 26¢° T°TS 8926" Z°L4 ot 9T 
LS¢ ° L°4g = 9S%° G* th COCs nee Che ett: 9°4G6 4h" 9°9S6 STh* L°2S ST GT 
9L¢° L°6e ¢l¢° 9°9h 4d¢e° C° th 994° 2°09 6S+° 6°6S Sor’ 2°Ss hT bT 
GOt4* Q°ch 404° G°0S 604° 2°84 = TOS* 4°S9 64° +. ° +9 T9t° G°9S Si cl 
get° T°9h 839 HE 4° 2°4S = LEH” 9°TS ¢¢S° G°69~—s-: T2S° 0°99 S§h° 9°T9 EL eT 
GLA 9°0S = 9zt° 9°65 =ORh°? = =6L 9G S° Z4°S<L = =O02S° 2°4d = =60LTS* 9° SY jee TT 
GeS° G°9S —S6SS° 6°99 045° 9°99 = zg" 6°TQ 919° €°08 9SS° 9°0d Or OT 
G6S° 6°29 969° G*tl 66S° Z2°OL 489° T°62  #©699° T°28 S0O9° 9°92 6 6 
439° Sir ee oe 9°Sg8 689° e°Tg 2S2° Or OGue cea. 6°46 T69° 9°28 Q ) 
66c° €*t3 TOS’ O*OOT. =TOQ i. -S* 46" -609" 4°SOT 409° 6°OT 103° 6° TOL Z ‘@ 
GO0°1 “9°GOL "COO Is 6 tel -COO;E  LiGll, COO'L vy Ol. 600 TL. 420C 0 ) 000: fk O° cet 9 9 
° ON 
9a/d d Ia/d d Ea d Paya d Ia/a d Ia/a d UOTIEIS aQny 
=O fe) a 8 fe) fe) fe) 

Ign°OOTt= 1 1,0°¢d= L A,0°6t+= I aeOnGh== al athe Ab eG 
etsd L°96=° 4 etsd /°¢tI= a etsd2: 2OT=d eted 2*2tt=q eTsds/tt= a etsd LeHtt= d 

are “ON 3994, II¢- *ON 389] _ “ON 399], II¢ “ON 3Sa] Ear “ON 38991 Pee “ON 4895 

he “CN edtsty Joy “Urq edt ceae 





| ee eS Ae 





A36 


Le 


GZ 9¢ €¢ cc 12 0c 


6l ST Lt ST cI 
adadnNN NOLLVLS 


al 


ned 
/ 
/ 
/ 


Ol # ISaL — — 
ll # LSad x 

6 # LSAL o 

8 # LSaL ° 


£ # LISaGL Vv 


JO TIONV ‘Z ALVId TuNsSsayd SII FTIZZON 


SaHONT xX 


0°t 


vA! El cl Il Ol 


o8l 2 = AONTOMAAIG 
“STIITAONd OLLVY FUNSSANd TVINANIWadKA ‘CZ ANNDIA 


0°T 


8a/a 





A37 


jee ren €°9¢ 
20% ° T°6¢ 
92c° 2°ch 
GG¢° 0°94 
6¢° O°TS 
O24° 4° 4S 
ehh’ 2°ds 
92° 4°T9 
88r° T°¢9 
TOS° 6° +9 
6TS° T°29 
245° 2-04 
T9S° 9°22 
9S ° 9°S2 
909° G°gd 
G¢9° 2°28 
999° c°98 
202° 6°06 
G42° 4°96 
692° o- cOu 
hhQ* 2°60T 
6T6° O°6TT 
000°T 4°62eT 
Sado “(Smid 
CUES 1s =°% 
etsd 2°*2 Tt =°g 
tr "ON UNnYy 


Coc G°92 gcc” 4° oh GT?° T° Te 
L62° 6°92 2S¢ ° 9°SH Lie L°¢e 
2ck” 4° TS gcc” 6°64 L9¢° €°9¢ 
GS¢ ° 9° #¢ Tt ° O° 4G 96° ¢°6¢ 
96¢° G°Qe O9+° 8°6S 475 +7° 6°?h 
6TH" "0H Lét* ee" e9 644" 4° th 
Gtt° ¢° oh 20S° G°S9 Th’ G°On 
Ogt* L°9+ ¢¢S° S°69 64° 8°Sh 
T6+° Q° lt cS ° 9°02 TOS° G°64 
TOS° Q°Sh 09S° ohare & OTS° 4° 0S 
22° 8°0S 49S° Oooo 9TS° O°TS 
G+4S° O°ES GQS° 0°92 THS° G°SG 
G9S° 0°SS 409° 4°92 6S$S° 2°Gs 
9Q9S° E225 ST9° 9°62 QdS° Teas 
€T9° 9°6S c¢9° T°2g ¢09° 9°6S 
G¢9° 6°T9 649° 2° HQ 0¢9° 2°29 
029° Z°S9 499° 4°98 649° T°S9 
GOd ° 9°89 902° 6°T6 969° 6°99 
94d ° 9°22 Ord° 4°96 Trd° 2°Eh 
€9d° e°9c 062° T° <oT 292° e2e 
THQ° 6°TR G4Q° 2 Ol THQ° T°€Q 
QT6° €°68 616° 6°6TT 616° 9°06 
000° T ¢°26 000°T T° OST 000°T 8°96 
8a/a (Sua Safa (BHn)d = Baa (BH) 
Ioft6L = d 1,9°SE- = Io9°de- = ab 
etsd 1°49 =°4 etsd ¢° TT =°g etsd 2°92 ae 
Tl oy ‘ON uny 16 “ON Uny oO -ON Uny 


GZ °onN etnStq JOl B3eq auNnsselg 


Q ATavVs 


90¢° G*Oc 
2s ° L°O¢ 
¢9¢° 9°S¢e 
6Q¢ ° 0°9¢ 
T2r° 0°6¢ 
Goh? €°On 
9S*° er ey 
¢Qh° Ashe 
C64" L°St 
¢0S° 9°94 
QTS° 6°24 
+S ° 4° OS 
29S° T°2S 
¢9S° O° 4S 
609° 7°96 
259° G°9S 
299° €°*T9 
COL. T°S9 
Ghd" 0°69 
292° oa tHe 
948° 7° 9d 
Teo ¢°SQ 
000°T 9°26 
Ca ae ua 

O 

I,9°0¢- = A 


etsd 2°0Q =" 


Il), "ON Uny 


*ON 


"ey 


“ON aqny 





YAGNNIN NOILVIS 


€l ‘21 # Sisal — — 
91 ¢ ISaL V 
ST # dSdL 0 
71 # LSaL ao 


ote’ % = AONASUTAIG 


JO FIONV ‘Z ALVId FUNSSaAdd ‘II FIZZON “STIIAONd OLLVa FANSSAad ‘IVINAWLUAdKA 69C MANOA 


SGHHONI - X 








A39 


09S° ="a/%a otyey eunssaig TeIOL g UOTIFIS pownssy . 


Ogg ° O° 2h o¢g° 0°6¢ bi Se O-2S ¢i2° 9°62 Ga2° B° co Ane G2 
rh 9°Sh Lat * G°Or G3 ° 2° 4S che” G6 Che CoG. ACZ hd 
QSz° 9°St 19¢° eer L9¢° 9°9S 93¢° tts LCs GC * 4g Acc C2 
GL¢ ° Q° hh Qd¢° a it €Q¢° T°6S G62° 9° 2e h62° T°S¢ AT? ce 
06¢° 9°64 46¢ ° T°9n 66¢ ° 9°19 QT ° 6° He Ci. vos HO? ie 
OTt* eS Ih? 4° Qt O2h° 6°49 Lec G° Ge O2¢ ° 2°6S HOT O02 
624° 9° 4S 9¢4° 6°0S cth* 2°39 Bes °9¢ Chg ° 6° Or “Ay l 6T 
ES; S°2S 094° L°sS 994° cet h9O¢° G°6¢ 99¢° 9°Ch ALT ST 
09° o°8S 99n° L°S 9Zh° G°ed cle ° €°On G2¢° Leh ea ZT 
Lo 0°09 TQ4° 2°9S 064° 9°S2 O06¢° C° eh 46 ° O-23; AQT JT 
T8t° 2° 9 064° G°hs ¢0S° O° GL OOn° 4° o + SOr° Leh 9T GT 
647° 6°29 805° 4° 6S hes: 2° OQ ten L°St ccth* #° OS HAST tT 
90S° 6°49 GeS° ¢*T9 Q¢S° T°S3 cot" 6°9h got" 0° 2S ST vi 
oe: 4°99 Qs’ 6°29 GSS° 9°S Qtt° 9° St TSth° B°ssS ATT All 
€SS° 4°99 €SS° 9°49 025° 6°28 GGt° 4° 64 6S4° Las 4T (BE 
OSS° G7 On 225° 6°99 465° 9°16 GOS° L° 4S 205° 39°09 AST Ot 
69S° O°ed €6S° 2°69 419° 6° 6 61TS° €°9S €2S° 4° 29 eT 6 
625° G°S2 219° yale 269° G°26 2G" 2°9S GuS° 0°S9 ALT -, 
209° eye 1¢9° gs 369° 4° TOT 69S° 9°19 Cus €*Q9 et re 
G29° T°Og 449° eo SC 929° G°*hOT 969° Lo aG 96S° meee ATT 9 
249° O°¢a €99° G*26 Cia: (Seo Ee T¢9° 6°39 4&9 ° a°G2 eal S 
¢49° €°OR 039° G°62 hyd 9°hHTT od9° 9°22 200. e7) HOT 4 
real i €° Th oy eas ota, H9¢° 5° oT Lice rae ¢1d° T°S8 OT ¢ 
928° €°SOT ¢¢32° 4°96 223° Suc | 213° T°9S Tie; 6°96 6 2 
000°T Cet O00" T 6°9TT P 000°T G*80OT 000°T €*é6Il Q if 
: ° ON 
Ra/a (2H )d 8a/q (44,)¢4 Ca/a (2Hu)d Ba/a (4H )d Raa (FHu)d uoTYFYS eqny 
1,0°66-= Ao2°on-= 1.0°26+2 i,921= 1o2°9nT= d 
etsd 2*2Ul= 4 eted 2°S0t= d etsd 2°GCt= d etsd 2°96= 4 etsd 2°90T= d 
or JON 1S 'o), Soil “ON 389] ON 399, et BON 2591 sar “O 4S4, 


a ee ee. eee 


92 °O" ain2@Ty aAOJ B3eC vINssadlg 


6 TIEVs 





Aho 


Ee cc re Od ot BT Lt QT ST TT 7 er TT OT 


YadWN NOTLVLS 


—e 
x Be 
¥ ' SS. 
SS 
} 7 oe 
& a 
ie ss 
a ¥ bene 
v 
a er 
% : \. 
Poe 
on «(CP . 
9 
OX 

~ 

NX 

aN 

a 
€l ‘2l # SISSL —— 


6T # LS&L Vv 
Ql # LSYL ° 
LT # LSGL a 


oce c = GONADHAATC 
dO WIONV ‘c GLWId ddNsSadd ‘II WIZZON ‘SH1IMOdd OILLVE SHNSSHdd ‘IVENAWIdddka °27 ddNold 


SGHONI - X 


O° O’d Ort 


d/d 





Aul 


492° 2°02 cee: Cake Aye Ge 
¢6e° Ete Qce° CE Gre’ T° th AS 2 4d 
9 Sa C°Se2 ¢G¢° L°o¢ T9¢* L°94 7.2¢ C2 
Oe ere OGs” e°G¢ gz¢ ° C° Qh Ale ce 
OE” 2°92 QQ¢° O°d¢ h6¢° 4° OS Ors te 
TOn* er ec Oh" 1562 GTh° 0°¢SS HOT Oc 
61” c° OS 624° 6°04 Goh? 9°SS AQT 61 
nt" Cee ¢Gt° eee 094° 3°86 ALT 8T 
cSt 9°2¢ 29° Lae 694° 6°6S ZT Gl 
G9n° G°eg 944° aoa cQh° OFS AIT val 
gL" C° He 98h" 4° On 264° 6°29 9T ot 
46t° 9°S¢ GOS° o° Qh OTS" T°S9 AGT HT 
real 6°9¢ #7S° 0°0S 925° 2°29 GT CT 
oe ae THS° oS THS° 1°69 ArT et 
HhS° ra ¢GS° 9°¢S GGG: Ore 4T U1 
Corer es)) 49° 6°0h = 2d" 9°4S gg Ene7 AE OT 
19S° lies T6S° #°9S 96S° 4°92 at 6 
06S° G* eh 709° 9°25 S19: Z°g2 ACT 9 
06S° G* eh T19° €°QG 049° bans @°. Zi Z 
¢09° G°oh Ceg" ¢°6S 76a O° 48 ATT 9 
QC9° 0°9 On9° Loa.S 489° oes, IME G 
c49° S°SH ¢29° 2°49 moc 9°68 AOT 4 
612° 6°TS Siz ¢°99 cee’ 4°26 OT ¢ 
G2Q° 4° 6S G2Q° Gee 929° g°SOT 6 e 
000°T saree 000°T 4°S6 000° T 6 cel 8 i 
"ON 


84/4 (7H. )d 8a/a (FHi)d 84/4 (2Hu)d uoTIEIS aqny 


2 0 - _O e _o 
TOT eye aL We cous AO°eh-= L 
16) oO 16) 
etsd 2°49O- d etsd ¢°HQ= d etsd Z°2Tl= d 
Iles "ON 4Sa] 115, “ON 3¢99L Iles “ON 389] Syren “ON 399] Il, "ON 3891 


-- ee oe a ee 


LZ *oN aanStq s0j ejeq aansseaig 


a 


OT ATAVL 





EC 


CC Te 


qa 


JO TIONV ‘Z ALVId FNNSSANd ‘SII FTIZZON 


0¢ 


4: 


61 


ST Ll 91 


qa 


€l°7@l# SISaL — — 
Cc# LSAL a 
lé# LSaL 2 


Oc# ISAL V 


cl 


7T tl cl Tl 


YqeHON NOTLVILS 


otz'% = FAONASUAAIG 


*“STIIAOUd OILVA FMNSSAad TWINININAdXA ‘QZ TANIA 


84a/d 





AL3 


neg ° L°6¢ GI¢° O°2¢e Oss° 4° SE Ane G2 
Cog® 1G; Gee” O°S¢ CH? Q°9¢ ACe hd 
6h ° O°o# L4¢° ¢° Go QG¢° 4° QE ACC v2 
G9¢° O°St 79g ° 9°9¢ He ° T°O# 412 ce 
Og¢° Q°OH 6Q¢ ° G° 6S 26S ° mes AO? Te 
O04" 6 °Gh Cth’ O°?h TT" Eat WOT Od 
rag te UG 9¢t,° Gh Lot 2°On AQT 6T 
Ont? 2°nS Q94° Q°Oh QGh° eden, ALT mai 
Cae Z°cs OZn° Q°dh 994° @ "06 ae reas 
9Sh° G°aG TR? 6°Qh T9t’ 9°TS AGT OT 
G9r° 4° 2S 264° 9°0S Ack 2°25 9T et 
O8h° T°6S 21S° 0°<¢s 10S 2° 4S AST tT 
26h" 4°09 O¢S° 6°SsS 62S° 2°9S en eT 
ZOS° 6°19 94S° SSS QhS° 2°95 “41 eT 
( 6 eTQeL 99S) TG T°s9 19S° Q°2S T9S° 6°09 4T wil 
TES° G°S9 €9S° 2°6S 99S° 9°29 oT en 
¢4S° 0°29 409° ¢°T9 909° 0°S9 oT 6 
9SS° 9°89 SHES SECeS AS Ss) AcT 8 
Sg As Leos Ces Q¢9° 4° $9 et Z 
TO9° 2° HG 59° €*° 49 G49° T° 69 ALE 9 
629° 9°92 249° 2°S9 999° €°od TT G 
G29° 2°SQ ¢29° 4° Q9 629: Cen HOT th 
9TZ° 4°93 9TZ° eral ole 9°92 OT ¢ 
929° 620i. + See? Q°Q 928° 9°88 6 Zz 
000°T €°oeT OOO°T 9° TOT 00C°T 2 208 g al 


“On 
8a/a4 (FB )d Sa/a (3H )d 84/4 (FH )d uotzeys aqny, 


10) ‘e) 10) 


d,S°é¢e-= & AS Hye L ta gies wb 
O O O 
etsd £°Q0T= d etsd 2°u6= d etsd 2°S6= d 
Bae i al . 
Go aon, Zio Ser ee Oe om tz con ysen Moz + on 3800 


Q2 °ON aInEty Joy Bzyeq aanssoig 
Tl SIdVL 








ALY 


7¢ 


ao of 


4 


o<] 


Ze 1c 0Z 61 ST | 9T CT val €1 ZI Il OT 6 
UAGWAN NOLL 
ie 
° oo 

<= 
V o™ 
° V Vv ere 

re] ee oe 
a o ° V ~~ 
a Vv ~ 
A Vvg Vv 
a @ ago VV eae 
Big tostaet a d/d 
a 
\ 
Lt # ISaL —— a 
o 
9% # LSAL © x 
CZ # LSAL Vv 
# \ 
92 # LSAL ° \ 
€Z # LS fa] 
\ 
x 
908°O = AONADUAAIG AO FIONV ‘2 ALVId \ 
AUNSSANd ‘II PIZZON ‘“SHIIAOYd OLLVA AUNSSANd ‘IVINAWINAdXA “67 FUNOIA \ 





ALS 


L0S° 4°09 
GTS° 9°19 
rcS°* 2°S9 
eeS° 9°S9 
O4nS° 9°49 
645° 9°S9 
T9S° LULZ 
945° 6°29 
€9S° 8°59 
665° Or Le 
n9G° oan 4 
609° 5°22 
G29° Ce 
649° ook 
G¢9° 0°92 
199° T° 6c 
929° (han Re, 
967° 6 °oR 
gee G°S2 
Te. G° 2g 
262° 9°06 
ver 6°eG 
GOxz #° 9G 
poo: 2G 
GOO oa So 
a/4 (30 )d 
Oo 
CLO T= 2 
0 


etsd L°S62 a 


il,s “ON 4SaL 


G2S° T° 99 
coG° T°59 
CHS° oO 
HhS° Cd, 
OSS ° eee 
6S¢° G° 22 
69S° O° 
HQ2G° CSG 
29S° O72 
G6S° eon 
009° Git 
TQ: 2°62 
G29° i be 
¢o9° i hares? 
G¢9° 4° e8 
959° €°SQ 
229° BRT e 
469° 0°06 
etc: ¢°26 
fen. 4°16 
GS" 0°26 
ree 4° TOT 
GOR° f° OT 
ore: OAT 
00071 CO 621 
Safa (FI) 
Oo 
aS eae dt 


ie) 
etsd ¢*TOl= d 


lee *“O'J jsoy, 


62 °ON aan®tq 10g e}yeq ainssarg 


é6S° 9° £6 645° 9°dg 
ods" —-h°66 645° 6°68 
18S" ars 265° 6°16 
2S" 0°96 965° 9°26 
46S° 2°66 609° G* Hb 
c09° Lae S| Biehs ees 
ATG ous eat HOO" 6°95 
O29") 26 SO) Seis) oor 
G¢e9° 9° 90T 659° 9° TOT 
SS eR Se ue Ce lO) 
orcs 2 COR UGor O°ZOT 
G49" O° QOT 6S9° 4° OT 
269° S*60T = 699" 6° €OT 
¢59° 2°60T 693° 62 COT 
059° T°6OT = 199" GF OT 
S205 6° TTT 329° 2°SOT 
n89°  gthTT) = $69" -9°2OT 
15 3°9TT 002° 2° SOT 
HTd° Soot CU LoOEL 
sec! eee eee top 
952° coer ote he cia 
eoae tan (ya O°T2T 
203° L° Het €03° Ley 
Pasd (3H.)d Bala (FH )d 
EGS eee 4 6°66-= 


ets L°tSte 


LE 


G2 


229° =a /8 gq OT}eN aANssadt THC], g UOTFEYO pounssy , 


nie 
ae 


"ON 38a] 


ol ATaVd 


etsd 2't2T=— 





"OY 48a] 


A 


Ne 


695° 4°9 
109° 4° 36 
419° 4° OOT 
419° h°OOT 
€29° O° 20T 
G29° €*20T 
6¢9° G*+OT 
259" Z-S00 
6S9° 6-20" 
199° Cee al 
799° 9° gO0T 
699° +> O0L 
139° Sie ae 
909° rere 
T9@9° Cr Ur 
909° T°HTT 
612° Caring 
veo olen! 
Ho 2° ee raal 
ced? 6°6TT 
092° 17 He & 
Te27 Gi Lek 
€09° Teer 
Paya M(H 
12) 
Oe co 


0 
etsd 2°22T= d 


ee 


4 


qsaq 


WH? Ge 
WSC tre 
Roe ra 
rte de 
HO? 5 Bes 
AGT O¢ 
ART 6T 
521 RT 
aml om 
HOT oT 
GT GT 
HST hT 
Gt eT 
AnT ean 
HT EE 
AFT Ot 
¢T 6 
ic T 5 
al c 
ATT 9 
TY G 
AOT 4 
OT 4 
6 re 
R eG 
"ON 
UOTRIEYS eqns 





Alit 


%¢ Ec CC T¢ O¢ 6l 


° pet een See 
° = w 
v 3 v $ ~ 
‘ $ 
° ° 
°o 
° 
oa oO a ie | 
aTINSSayd SII FIZZON 
O°€ 


ST LT 91 SI va el cl Il Ol 
YHaWAN NOLLV.LS 


ie 

oo 

~ 
~, 
° = 
: o 0 & 8 oe 8 aS 

= o o =+G Soles ae 

o “XN 

NX 
‘XN 
L@ # ISaL —-— aN 
\ 
Tl # LSA °o N\ 
BS 

oO¢ # ISAL Vg \ 


60 # LSal ° 
8c # LSaL a 


908°O = AONADYHAIG dO WIONV ‘Z ALVId 
“SHIIAOUd OLLVY SUNSSHYd ‘IVINAWINAdXaA SOE RINDIA 


8a/d 





ALT 


(2T 8TQeL 295) 


42S" 4°49 
TES° €°S9 
OnS° S°99 
Czy 4°99 
QS" G°29 
9SS° Se els 
99S° 6°69 
695° Sete 
295° ced, 
G6S° e°o' 
¢09° ine 
TTS* Ga 
Bc9* eas 
G¢9° ter 
9c9° G29 
6S9° Le 
929° 2°S8 
£69" +°S8 
Ce (re 
6c2° J" 65 
962° T°S6 
¢6h° ¢°96 
GOR" -° 66 
828" Z2°Qot 


00c’*T ee 


8a/g (48. )A 
0 

AS Te-= dL 
@) 

etsd C *O6= a 


Tlre -oy yan 


Lc9° 
9¢S° 6°92 
hhS° 7% Oe 
255° H° TS 
2SS° 7°TQ 
09S° G° 23 
G9S° c°S8 
Sia. 8°48 
88s" O98 
T6S° T°2g 
96S° 6°29 
09° 2°99 
tee T° 06 
969" eb 
S¢9° G°6 
829° C°n6 
655° 0°26 
8ZS° 6°66 
469° ree rash: 
rl 0° SOT 


Bea” 277 OW 
T92° © Zit 
622° 6° nT 
008° OC OLE 
C28" 0°62T 


84/4 (28 )d 
oO 
d,0°6¢-= L 


0 
etsd ¢2°Stt=d 
Tlo¢ “ON 3S9L 


ee 


OF “ON ounety soy BYEq suNssaid 


€T TIGYL 


9SC° Coe 
02S° G°¢cd 
(Seles eae 
69S° 6°S2 
409° eric 
¢19° 0°62 
Q&9° Vibe 
949° C°S 
259° u* He 
¢S9° T° +Q 
0S9° 9° SQ 
649° LER 
759° ¢°H8 
TS9° 6° EQ 
hho? jose 
799" Q°SO 
7g9° Cae 
869° 0°06 
Cis 2°26 
AE 4° 46 
QG2Z° £26 
G27% E-OO1 
405° oc OL 
Sag (3H.,)d 
i) 
x) a aC 
etsd 2°10T= 4 
Tins +on ased 


9T9° 
9¢9° 
949° 
8S9° 
659° 
959° 
+99° 
229° 
699° 
249° 
999° 
Coos 
999° 
Q9°° 
¢S9° 
699° 
S89° 
G6Q9° 
GTd° 
6¢2° 
294° 
GeL. 
GOd* 


8/2 


d 


© get is 
az 
G°9/, 
2°92 
8°94 
7°92 
CG 
¢° Qu 
2 Oe 
$ ° QA, 
2°22 
oce 
QrZG 
6°92 
0°92 
G° dd 
ZL. 6c 
Docs 
¢°oR 
6° 4Q 
Lk Oo 
9°06 
76 


(2H )d 


"6G-= 1 


etsd L*6= 4 
Tl oz: on 4895 


witiC 
AS e 
Kec 
¥.1¢ 
Ora 
Ao 
AST 
nal 

LT 
n9OT 

ST 
AGT 

GT 
AnT 

vale 
HS T 

eT 
ACT 


UoT}EYS 


°a/2a OT}ey aInssadg TROL g uoTt4AS Founssy , 


OT 


a Ae = UN «OD C= 20) O* 


Orr 
eqny 





140 


P - ATMOSPHERES 


FIGURE 31, CONDENSATION DATA. 
LOCAL CONDITIONS AT ONSET 
OF CONDENSATION. 
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FIGURE 36, ONSET OF CONDENSATION, NOZZLE II, 


COMPARISON OF THEORY AND EXPERIMENT. 
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APPENDIX B 
REVIEW OF NUCLEATION RATE EQUATION DEVELOPMENT 
The following is a review of the classical nucleation rate equation as 


developed by Jd. Frenkel?) , 


Although Frenkel was not the first to investi- 
gate nucleation theory, his result is not greatly different from earlier 
authors mentioned in the text, and his method provides relatively straight- 
forward consideration of the important aspects and assumptions in classical 
nucleation theory 

Frenkel's droplet model was based on Volnen eee model of drop- 
lets as generalized molecules, neglecting intermolecular forces, and his 
method of development was based on that of Pevaowmien co He showed that 
the distribution of clusters of vapor molecules, resulting from the kinetic 


interaction of the molecules could be described by a Boltzmann distribution. 


That is, that (in equilibrium) 


~A/KT 
N(g) = Ny exp (A#/KT) (B-1) 
where Be = total number of molecules 
®@ =H ~- TS 


Ao Sa >, = change in free energy of a group of (g) free mole- 
cules when forming a liquid embryo. 

Several simplifying assumptions of debatable validity are made at this 
point. First, the droplet model is taken to be a liquid sphere, having sur- 
face energy described by the flat-film liquid surface tension (or some cor- 
rected function thereof). Second, the total change in free energy is assumed 
to be the work of formation of the surface of the droplet, plus the change 
in free energy resulting from the change in state of g molecules condensing 
from vapor to liquid, with the spherical droplet assumed to be stationary 


in the vapor system, and having no rotational or vibrational energy. This 


is the important distinction between classical and more recent "gassified" 





be. 


droplet models of nucleation theory. Third, the droplet is assumed to form 


(32) 


at environmental temperature, an assumption which may also be invalid 
Using these assumptions, and defining $ as the free energy per molecule, 


Ad can then be writeen: 


Ad ao, = $,,) + nro (B-2) 


~ T(AS) + dnr2o 


q 


= ek on ae lirr2o (B-3) 


ao 


where p. equals the flat surface equilibrium pressure of a formed embryo of 


(g) molecules and p is the existing local pressure. 


Hence, 
(Canes eran cing os aes (Beh) 
E O - i KT = 
paying that: 
eae. 
gm = amr°py (B-5) 
yields 
= RB _ do (3gmy2/3 
N(g) N exp { g Xn >. ~ KT hea) } (B-6) 
2 Dea e ees 
Ny exp { g &n D. Co) g } 
where a no (BE)23 


Inspection of Eq. (B-6) shows that a reasonable equilibrium distrib- 
ution results when p is less than p> that is, in a superheated vapor. How- 
ever, in a supersaturated vapor, an unreasonable and, in fact, impossible 
equilibrium distribution is predicted as seen in Fig. (B-1), since it pre- 
dicts N(g) >» as g>o . Acknowledging that an equilibrium distribution 
of drop sizes is unreasonable under supersaturated conditions, early authors 
treated the system as one in a state of quasi-equilibrium, where droplets 
are imagined to be removed after reaching some size larger than the critical 
radius, r*, which occurs at the maximum of the supersaturated free energy 


change. This maximum in A (minimum in N(g)) occurs when; 
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Bu. 


aa eS 
(Ad) =O -—-KT mh 2 4 5 Ef, / 
ty Pe 
| 2,-/3 
Kv! &n & = er 
Dee ie 
dee aa 
Ad = —-€ — Tr“o 
max 5 : > 2 
Tnen 
| yy ”) 
Ad = — okT & Se + lt 25 = = 0 0 ae 
max eae : a (B-7) 
Hence 
3 Pp One? 
(= ae =— T7 0 
KT(3 pes ) &n ieee r 
mo 
ee) Sune greta dea 
ine Omer (B-8) 


Equation (B-8) is the Kelvin-Helmholtz expression describing the conditions 
under which a liquid drop may remain in equilibrium with its vapor phase. 
The radius (r*) and the number (g*) of atoms of a critical embryo are thus 
determined. A variation of r (or g) in either direction results in a de- 
crease of free energy from the maximum, and consequently the drop is in 
unstable equilibrium at this point. 

In considering the quasi-equilibrium state, and ignoring interaction 
between nuclei, it is said that nucleation rate is the net rate at which 
particles of a given size (g) are formed. T>% maintain quasi-equilibriun, 
it is further assumed the net rate at which (g+l)-sized droplets are formed 
from (g)-sized droplets is equal to the rate at which(g}sized droplets are 
formed from (g-1)-sized droplets, and so on. That is, nucleation rate is 
independent of g, in order that the number of drovlets of any size remains 
constant and the rate at which droplets of larger size < (Gaeined <2 be 
removed is also equal to the nucleation rate. Correspondingly, it is im- 
agined that single molecules are added to the system at a rate equal to the 
depletion due to nucleation. 

Using this model then, and assuming that growth of nuclei is due only 


to interaction with free molecules, the nucleation rate is 





i 
J(g +1, t) = B s(g)x nlg, t) - ylg + 1) s(~e + 1) n(g +1, t) 


(B-9 ) 


where 8B = rate at which molecules strike unit surface area. 


y(g) = rate at which molecules leave the surface per unit area 
n(g,t) = distribution of non-equilibrium drops 
s(z) = surface area of a g-sized drop. 


Assuming that no effects of droplet curvature exist, and neglecting the 

cluster motion, 
8 = Canc Py IZ 

Also, it is said that if nucleation rate were zero, then the non- 
equilibrium (or quasi-equ#librium) distribution would in that case have to 
be equal to the equilibrium distribution N(g). This is easy enough to see 
in the superheated case, but perhaps a bit difficult to accept when N(g) is 
the fictitious mathematical result corresponding to supersaturation. 


Applying this, when J(g, t) = 0 


Ye) set) pve = Bese) “nte) 
when n(g) = N(g) 
Hence, 
y(g + 1) s(g + 1) N(g + 1) = B s(g) N(g) 


and Equation (B-9) becomes 
= (sre eal woe ty) ~_ 9ry 
Ig, t) = 8 s(g) N(e) ( 7S ented (B10) 
Here some rather notable simplifications and manipulations facilitate 
a solution to (B-10). For large values of (g) where these functions may be 
treated as continuous, this is approximately: 


J(g, t) = 8 s(g) Ng) = (2) 


dg N 
Hence, since it is assumed that J(g, t) = J(t) 
Oe ne eee Beal 
} My) = -&F Sle) Ne) ae 


Using the boundary conditions 










_ 
mio an heoe bere olde, Me ae 

- | >= 
2 








_ _ 
peta ("=v 72 ce ae - ~M 
re! 
—_ 
7 Mill “The bitin taep 7 } 
a 





7 - 
7 AS 








i in) pre 


- 
ie : rhe ud 
_ 
ar is aid sit ‘ 
7 
\ i = . 
ee 
_ 


nigl 


B6, 


(yO ete At) Lin 7 = 0 
ace 
One obtains 
J = + (B-12) 
fat sN 
Since ee ae xp (= = | has a sharp inaximum in the vicinity of g = g*, it 
N No” KT zs 
is possible to approximate (B-12) as: 
x 
- Beste") N, 
fyexp (A6/KT) dg (B-13) 


Using Equations (B-5), (B-6) and (B-8), the exponent may be written: 


Ad p 273 
pall a + 
OS ae 


2mo ing p 3m_ 72/3 ee) 
seu pkrr® J * kr inp J 


A [ 3(e/e*)2/3 - sale 





where 
4, = om aye 42/3 ree 
pei sane 3KT 
r*¥) N 
J = ix exp {A [ 3(/a*)?/3 - 2(e/e*)] } ae (B-1)) 


This integral occurs in the Becker and Doring, and Frenkel derivations, 

and can be solved approximately with the introduction of 
9 = (g/g*)1/3. 1 (B-15) 
Since the integrand of (B-14) has a very sharp maximum an at g = g*, the 
important values of 2 are small, and the limits of the integral may be taken 
at +~ and -©, 
Kiso. ‘since (g/e%)2/3 = 92 + 22 + 1 
= 3(1 + 2)2 g*¥an 

And 3(g/e*)2/3 - 2(g/g*) = - 322 - 224+ 1 

The integral may be written approximately 
3e%e f° (1 + 2)% exp (- 3An? - 2A03) an 


= 3p%e" J", exp (- 322) dg 





By. 


: 2 
= 3g%eM(n/3a) / (B-16) 
Hence, the nucleation rate finally becomes 
& 
_ Bs(g_)Nop A_ 53/2 
= ei al exp(-A) (B-17) 
with 
No = p/KT 
8 = p/(2mmkT)!/2 
*)\2/3 
Ei pene alae 


and with the definitions of e« and A, (B17) can be shown equivalent to 


#9 
= ee 20,\1/2 ~inor 
g= (B) 2 (2)! /2exp(i—-) 
KT" p 3KT (B-18) 


This is an expression of the number of nuclei formed per unit volume 
per unit time, due to the existence of a non-equilibrium distribution n(g). 
By assuming (J) was independent of (g) as Becker and Doring did, it was 
evaluated using only the boundary conditions of n(g) and the properties of 
Ad(g). 

Discussion and criticism of Eq. (B-17), and of the importance of its 


various constituent values is contained in the text. 
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APPENDIX C 
REVIEW OF DEVELOPMENT OF DROP GROWTH EQUATIONS 


The following development follows the method of Oawatitcn 2. 


in consid- 
ering the necessary conservation of mass and energy of growing droplets. 
Included are a correction of average incident energy to be ORT rather than 
3/ORT, and with generalized non-unity values of both mass and thermal accom- 
modation. An important assumption is that molecular velocity distribution 
are Maxwellian before and after collision and accommodation with the drop- 
let surface. 


Consider a drop of temperature T, in an atmosphere p and T. It has an 


D 


incident mass flux equal to 


(2mRT)1/2 


with a proportion & (condensation coefficient) condensing and the rest being 
reflected. Since the drop is very small, its temperature is assumed to be 
constant throughout. The mass transfer from the drop is equal to that which 


would be incident in an environment Ty and Py» the equilibrium values of 


temperature and pressure, where 


omo 


Qn 
Peat Tp) 2 eae 


Hence, the mass flux per unit area passing out through the surface is 


EPy 
Thus, the conservation of mass may be written 
dr E SPp 
po, nr? cae) = tues Eoeesar yea Cart )177 | 
or ar E Pp 
ae “(acne | TZ 7 172) aly 


assuming that & is constant, or at most, a function of the drop temperature 
only. 


Considerations of energy conservation are made after first defining the 





G2. 
the thermal accommodation coefficient, involving the change in temperature 
of reflected incident molecules. 


Gus ‘vrefl. re (Thermal Accommodation 


C-2 
ue 7 Coefficient) ( 


The energy flux to the drop is, per unit area 
Mr Aaa se Che 
(areyt7z [ eRe | 
ana the energy leaving, upon reflection 
(1 - &) 
—= + — 
(2nRT) 72 {eR [ T + a(Z) 7) 
In adaiiiton, the evaporating mass flux carries energy 
Ep 
D 
(2nRT_) 1/2 l a ] 
D 
assuming the evaporating mass flux also to be Maxwellian. 


The rate of change of droplet internal energy is 


aT 


ar 
rae 3 ae eet 
cP kame: Har’ Fe 07, Yen 
£ ee al a 
Ven Dep y-l at ee 


Hence, conservation of energy requires, 


\ sa dr 
Lien eet) Fe ee Onaga Cease ee 7 
3 rpc aE hor ab or fo —y Bt he) 


g(2RT) - £8, (2RT,) - (1 - &) gar [ T+ off, - T)] 


ies) 
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or ro,c dT, boon T 
Ed oe ee ee la Seo) 





p 
TDS SOR ean Gina 
7 Rr ac [ vero ee Neg ] 
Pp 
Dire eee D 
(iderds B= (oaRT) By = (oapt_)1/72 


D 
An order of magnitude analysis shows that the left hand term is small 
compared to the last term on the right hand side. Taking typical values; 
r~10 7 em fe 6 Te 3 
fg Pp 
c ~ .3 cal/em °K R ~.045 cal/gm °K T — 250° 


yy - 1)~h dT, /dr ~ 7 x 10) teem 





eS. 


re oa 
a as a = .008 ~ 10-2 
fr 
(3RT)(~--—) (== - 1) 
° ee a 


this signifies that the energy required to change the droplet temper- 
ature during the expansion is much less than the latent heat imparted from 


condensing incident molecules. Neglecting this term, the energy equation 


becomes: 
Baal ab 0 
a eee ae Z ae das OS aly eee es 
O= cel - o> wp) - atl - €) ale - 1) +B a oe 
"te 
where N= CT 
p 
Since 
ieee 
as 0 (B - 8) (cay 


The energy equation can be further reduced to 


sianaih dl B 
02 (1-2) - oft - 1-1) + $0 - SDA - 


Core) 


Solution to the drop growth equations involves an iterative determin- 


ation of drop temperature, T_, from Eq. (C-3) and the resulting value of 


De 
growth rate from Eq. (C-1). 
An interesting indication of the importance of considering both mass 


and energy balance may be seen by closer inspection of Eqs. (C-1) and (C-3). 


It is seen that the requirement for a drop to grow is that 


B B 
Chloe ae) pee) Ro 7 
at * 9, (1 5 oe ime } (c-h) 
Re-arranging Eq. (C-3) into the form: 
Dp Py kr 1] <2 - ob -2y2- ye HY A- 
B iq 2 eal 7 7 E 7 ie ae 
we ohtain 
T 
eS peas Oe 1 ie alee 
By ee ee ee oe) ae 
Be T 


>, 


/ 





Oue 


It is clear that whenever T, > T the denominator of Fa. (C-5) is always 
creater than the numerator. Hence from Eq. (C-) the droplet will grow. 
Even in' the most unfavorable case, that of a = 0, where the dron tenverature 
becomes much hotter than the surrounding environment, the value of 8/8 will 
remain less than unity and the drop will grow. This is really a nroof of 
the intuitive fact that a droplet of size greater then critics3l radius must 
continue to grow in a supersaturated environment. Conversely when TT, 


6/8 is greater than unity and the drop evaporates, sunnorting our knowledge 


hae 


that evaporating liquids are cooler than their superheated envirenmtents. 
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APPENDIX D. 


MODIFICATION OF FOREGOING GAS DYNAMICS, NUCLEATION AND DROP GROWTH 
EQUATIONS FOR DEPARTURE FROM PERFECT GAS DURING AN ISENTROPIC 


EXPANSION. 

The following corrections introduced into the perfect gas relations to 
account for departure from perfect gas behavior are based on the empirical 
observation that compressibility factor z = Pv/RT is almost constant during 
an isentropic expansion. This is illustrated by Table 15 which lists the 
values of P, T and z for all of the condensation tests of this investigation. 
Changes in z generally encountered were less than .5% from stagnation press- 
ure to maximum supersaturation prior to condensation. This in effect allows 
introduction of "effective' perfect gas relations, 

Py = R' T where R' = 2R (D-1) 


and the gas dynamics are affected accordingly: 


Tds = QO = ae - vdp 
fe pdv . vdp _ 
p RZ + 6 vdp 
= aY| Py Sp} Yet op 
Y-1 p p Y p 
Hence = dp | ap | =) _ ap 
@) e p 1 z =) ee D (F1) (D 2) 


Hence, the isentropic relations of a perfect gas may be modified by replacing 
wih y yas defaned in equation (D=2). 

Similarly, the free energy of formation of a nucleus is modified as 
follows. The free energy change from vapor to liquid of (g) molecules (added 
to the free energy of surface formation in the classical droplet model) is 

dd = dH - d(TS) 


If it is assumed that the nucleus forrs at environmental temperature. 
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-TdS = vdp = 2RT- - 





D2; 


go, 


oo - 6, = zgkT &n — = ~ zgkT fn 2 
Py Poo 


or, the total change in free energy of formation of a classical model droplet 


becomes 


95 =— Oy = =ZekT cae + tr? o 


oo 


By proceeding as in Appendix B, it can be further shown that 


and 


(ee eee 
De P. 2KTr* 


These, plus the corrections to the drop growth equations, may be 


introduced to the perfect gas equations by the simple substitutions 


eee 
eer (D-3) 


' ny 
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In a computation involving condensation of a "real" gas, it is necessary 
to first determine the real isentrope, determined through use of the equation 
of state as mentioned in the text. 

From this, by use of the equation of state and Equation D2, respectively, 
z and Y' are established at all points of the expansion and it is then a 
simple matter to determine the value of y to be entered into the computation. 
The value of z is taken as the mean of stagnation and condensation values 
(indicated in Table 14) and a linear, temperature dependent approximation to 
Y between the same two points is made, 

Y= Y¥. -y¥, (°C) 
a b 


where values of us and Y, are given in Table 15. 


b 
Consequently, the computation scheme treats a value of y which will yield 


the real isentropic exponent y' when treated with the corrections (D-3) and 


(D-4). This matches the isentrope prior to condensation with the best 





Ds 
estimate of local conditions and calculates corrected nucleation and drop 
growth rates within the theoretical framework of perfect gas. 

Perfect gas calculations may be made by simply setting z = 1.00, and 
entering the desired value of y= Y_, Nes oe 

It might be pointed out that the compressibility factor begins to vary 
more rapidly after condensation commences, and that the computation scheme 
ignores this. Some error in profile shapes may result, but probably are no 
more in doubt than those resulting from perfect gas treatment. The value 
of this treatment lies in its attempt at establishment of the best estimate 
of conditions and behavior at the point of onset of condensation. 

Another possibility for the consideration of real gas behavior in 
condensation would be to completely rewrite the entire processes of nucleation, 
drop growth and gas dynamics in terms of an equation of state, as was done 
for the non-condensing isentropes and total pressure computations. This 
could only be done by greatly increasing the complexity of an already complex 
computation and at the cost of excluding the many perfect gas relationships 
which so greatly simplify the mathematical treatment of condensation in gas 


flow processes. 
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APPENDIX E. A CLOSED FORM SOLUTION TO THE MOMENTUM INTEGRAL BOUNDARY LAYER 
EQUATION INCLUDING APPROXIMATE EFFECTS OF STREAMWISE DENSITY 
VARIATION 


In most cases involving the study of boundary layers, effects of changing 
density are negligible. However, in such cases as nozzle flow, where sonic 
velocities are approached or exceeded, acceleration forces due to the presence 
of the density gradient become appreciable, and near the throat of a supersonic 
nozzle are of the same order as those due to area change in their influence 
on main stream velocity and hence boundary layer growth. A study of the 
momentum integral equation has revealed that a closed form solution is 
possible for both laminar and turbulent flow. In the laminar case, an 
extremely simple form of solution is possible, enabling one to determine 
estimates of boundary layer growth and depression in supersonic nozzles, 
using only slide-rule calculations. 

If it is assumed that density and pressure vary in the direction of flow, 
but not across the boundary layer, the momentum integral equation may be 


written as follows: 





d : dU : 
ax ? i (i u)udy + fevame (U = u)dy = TV=0 (E-1) 

O O 

NE 
—_——+ => 
x 
For the general case, it is assumed that 
Uy i 6) (E-2a ) 
1 

2 /n 

oe — ee (=) (E-2b) 


Then, Equation (1) may be reduced to: 
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d 2 aU Vv vee 
CT a 64 aU. 
DF (osu ) +d, edu = a, (=e) 


our (n=2)) 


where bj, by and a, are dependent upon the assumed profile. 


Equation (E-3) may be manipulated to appear in the form: 
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(Oy a Namie ee mes 
where ; D1 . a 
Equation (E-4) may be put into integrable form by using the integration 
factor: 
fPdx 
mG = shag = ce So) oe 
1 U fe) 


m= U? p2 (E-5) 


or, more conveniently m = (= ) ca 
+) 


and the new form of the differential equation is: 


a (m2) = ma saad 
il 
Ae 
where Q= C, (ais? 
Hence 
] 
d ica en 
5, més“) = C, m(=) Ag teeLe 
n-1 l n-1 
= ©, eee 2n (2y/n (m6“) On 
-1 
== n+l ] 
and: (m*) ep f(ms2) = Com en - ve 
dx U 
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Integrating, one finally obtains: 


nt | x 
a= C (n+1) n+1 j 
C 2 — Vv 
(m62) <" x, Qn Ix m en ae dx (E-7) 


Equation (E+-7) is the most general form of solution to the momentum integral 


equation, approximately valid for a compressible gas in either a laminar or 


turbulent boundary layer. 


It should be noted that the integrating factor m can often be more 


conveniently expressed as follows: 
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On amore conveniently: 


A 0 ] 
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mo 5 es na (E-8) 


Laminar Boundary Layer 


Considering Equation (E-7) for the case of a laminar boundary layer, 





where “= J, we immediately have: 
x 


x 
jms? | = mn m(77) ax (E-7a) 
O O 


Assuming that viscosity U is nearly constant and substituting for m with 
Equation (8) yields: 
Cu x m 
62 = 20 _ + —2 62 
m J moa) x m 
O 
A 0 A wie 1 6 
ee ce ~ ) —-) dx + — 62 
uw (F,) ¢ (=) JourGee acer eos 


Further, for a nozzle: mass flow rate = constant = PUA = P*U*A* 


Hence; 2 
Cou os oi =e oe 1-C, C=C) m 6 
ogee ee (nn gas (BN agg eo Rina) 
pee *A* ) 2 p m 
0 x O 





E4. 
Equation (E-9) is the specialized form of Equation (E-7) for a laminar 
boundary layer. Reviewing Equations (E-3) and (E-4), it is seen that the 
constants C, and Cy will depend on the polynomial used to approximate 
the shape of the boundary layer velocity profile. 
For flow in a nozzle, (A/A*) is known and C1, can be easily 
determined. Hence, the function: 


1-C) 2=-C) 


p 
Gad) (E-10) 
O 


may be also easily determined and plotted. The integral may be determined 


by the area under the curve F, and at any point, x: 


x 
ae co ee Oe ne 
Xx m 0 


6 = K m i F dx (E-9a) 


Pa 
where KR, = pey* 


Equation (E-6) was applied to a nozzle of known geometry, using the 


assumption: 


u 2 (x) 
eg ag 
9 


C6 


2 


for which (C, = 
At first glance, the assumption of a laminar boundary layer may seem 


(R1) has shown that turbulent boundary layers 


unjustified. However, Launder 
under sufficiently negative pressure gradients will revert to a definite 
laminar form. Also, for the conditions of this calculation the observed 
boundary layer at the nozzle throat was such that the boundary layer Reynolds 
Number was about 2000 at the throat. This is less than the transition value 
and would indicate that a reversion to a laminar form could have taken 


laces Hence, the laminar assumption is not unreasonable. 





ape 
Calculations were made according to Table I for CO, at stagnation cond- 


itions of 84 psia and Blowin according to the corresponding values. 


-8 a7 
F (A/A*) ~ (p/p) p= 8h psia 
30" 


O 
ee =10.8 


Ky for CO> at T 5 


Results of calculations are shown in Figure El. 


At the conditions noted above and for the nozzle in question, limited 
interferometer data provided visual means of measuring the boundary layer 
thickness, and four values are shown in Figure El. These observations are 
estimated to be accurate +15%. 

Also shown in Figure 1, to assist in estimating the relative size of 
the boundary layer, are those values of uniform flow flat plate laminar 
boundary lavers calculated for the densities and velocities at the entrance 
and: throat of the nozzle, respectively. 

It is seen that the limited experimental data seems to agree well with 
the approximate calculation of boundary layer thickness. More experimental 
data is necessary for a more thorough evaluation. 

The extremely small boundary layers resulting in the pressure gradients 
which occur when compressibility becomes important make it appear that 
further interferometer observations would be an advisable method of gaining 
such data if it is not readily available in the literature. 

The fact that the past "history' of the boundary layer has a greatly 
decreased effect on the value of boundary layer thickness in the region of 
rapid acceleration is quite apparent both from the curve in Figure 1] and 
from inspection of the equations. Equation (f=9) shows that the value of 
62, which exists at the beginning of integration must be multiplied by the 
Frauvo mo/m Hence, if the value of m increases appreciably, as it does 


when flow area A and density p decrease, then the term (m_/m,.) 5% contributes 
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little when added to the integral term. 

The fact that one is limited to an assumed boundary layer velocity 
profile during an integration does not restrict the possibility of breaking 
the flow into regions, and applying Equation (E-9) to different velocity 
profiles (i.e. different constants C; and Cy) over the various intervals of x, 
with relative ease. 

Turbulent Boundary Layer 

It is apparent that Equation (E-7) is readily applicable to turbulent 
boundary layers in its present form, merely by selecting the expressions 
desired for Equations (F-2a) and (E-2b), that is, the velocity profile and 
wall shear dependence. 

However, Since most persons are accustomed to referring to turbulent 
boundary layer behavior in terms of momentum thickness 6, and shape factor 


H = 6*/@, it is worthwhile to re-derive Fquation (E-7) in terms of 6 and H. 
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Hence, it is seen that: 
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Integrating yields: 
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Inspection will also show that Equation (7b) could have been obtained 
directly from Equation (7), using the relationship between the two inte- 
grating factors, and between 6 and 6. Note that C, = LG; (ea) he 


Substituting for m in Equation (7b) yields: 





nt x H+2 n+l n+1 
H+2,..\———= neti en oU —= H+2. .—— 
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where n= (H + p)dn_+ 1) a i 
Nn Nn 


n+1 


C= Gir a 
O 
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O 
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O O 
For the case of negligible density gradient Equation (E12) reduces 


tor nquat tone 22.0) "in Schlichting *?) , 
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for which 
Vy = 820 ft/sec. 


jee 325 lom/ft> 
Rote. 67. ft 


Fence, the estimated pressure gradient over the top ef the boundary layer 


at the throat is 


2 : 
(Cl oy a a oa? lpi psi 
an PR 1900 43 Jose as 


The side walls are approximately 1" apart. At a point midway between 


them, the pressure gradient is zero. Assuming the pressure gradient is 


linear, 
ig +h, 33 
az; 
cP = 33 + 8.662 
-—, 33 
z= 0 zZ = 1,0" 
And: 
ow Py pe PY 7 4.332 + 4.332? 
en em nese UnSEIe 
ne Py 2 3 O 
Sete il 
Hence: eG De = a psa 


This is of sufficient magnitude to offer a possible explanation of the 
less than 1 psi difference between average mainstream and measured wall 
pressure at the throat. 

It is noted that greater rates of curvature of the boundary layer occur 
upstream of the throat at relatively negligible velocities. Downstream, 
where velocities continue to increase the boundary layer has very little 


curvature, 


ls 
APPENDIX G 


DISCUSSION OF INTERFEROMETRY APPLICATION TO TWO-DIMENSIONAL 
NOZZLE FLOW 


Extensive discussions of the theory behind interferometry techniques 
are numerous. References (14) and (36) are included for the reader who may 
desire additional detail. 

Basically, an interferometer consists of a monochromatic source of 
light which is focused into a beam of nearly parallel rays. The beam is 
then split into two beams traveling separate paths of identical distance 
after which the beams are rejoined and focused on a viewing, or photographic, 
plate. An interference pattern of parallel light and dark bands, as seen in 
Fig. 12 is formed when the two beams are rejoined due to the fact that the 
various rays of each beam travel slightly different distances depending on 
the geometry of the mirror arrangement. 

One beam is passed perpendicularly through the test section. To com- 
pensate for the change in optical path length caused by the glass side walls, 
"compensating glasses of exactly the same thickness, must be placed in the 
path of the other beam. During testing, the only change which occurs is the 
density of the vapor in the test section which influences the velocity of 


light passing through it according to the relation, for two-dimensional flow: 


SX 


re ete it 
where n is the index of refraction of the gas in the test section, 2 is the 
width of test section through which light passes, A is the wave length of 
light being used, and S is the number of fringe shifts which occur at the 
point in question due to a change in the index from ny ce n,- 
The index of refraction is related to the density by the relationship 


Go =n-1 


where G is the Gladstone-Dale constant, fixed for any gas. 
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The technique used in this investigation was to first photograph the 
interference fringe pattern for no flow and then take another photograph under 
test conditions. The fringes of each were then numbered from a common point 
which effectively represented stagnation conditions. For each distance x 
it was then a simple, though time consuming, matter to determine the fringe 
shift S and thence to calculate the density ratio. Photos used in calcul- 
ating the density curves in Fig. 15(a), (b) and Fig. 16(a), (b) are shown 
in Figs. 11 and 12, respectively. A graphical illustration of the determin- 
ation of S for test 1 (Fig. 15) is given in Fig. 13. Test 2 was done more 
precisely by measuring the distances to each fringein the no-flow picture 
as measurements in the flow picture proceeded, since it was detected that 
slight departure from perfectly even spacing of the no-flow fringes occurs, 
probably due to minute differences in thickness of the optical glasses. 
(Note differences in Tables 3 and 4.) 

Although interferometric methods did not succeed in producing measured 
effects of condensation, as mentioned in the text, the results obtained for 
non-condensing flow does offer a recommendation for possible further appl- 
cation in this direction. One important unfavorable aspect of it, however, 
is the extreme amount of time and effort required to reduce the photographs 
to meaningful data. Each of the two tests presented here require alignment 
of the interferometer with the nozzle (a process which may take as long as 
two days, after which, however, it is permanent), photography, development 
of negatives and prints (another day), measurements from prints, and calcul- 
ations based on the measurements. If pressure measurements are taken, coor- 
dinating the location of pressure and density measurements also requires 
painstaking accuracy. 

Such a laborious process in reducing the data is still predicated upon 


a two-dimensional nozzle and becomes prohibitively difficult for three- 
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G3. 
dimensional flow. Reference 36 involves the analysis of one three-dimen- 


sional nozzle photograph. 
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APPENDIX H 


COMPUTATIONAL PROGRAM 


Comparison of Equations of State 
Isentrope and Enthalpy Calculations 
Normal Shock Calculations 
Nucleation Contours Calculations 


Condensation Calculations 
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S18 iat GO ea el 
IF (e¢9999*VP-V) 159265926 
IF AK=2) 27328528 
K=2 
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LABEL 

SUBROUTINE PRINT (NeXOUTsYeDY9J) 
PRINT FOR ISENTROPE CALCULATION 
DIMENSION XOUTI 200) sY¥(1l)s0Yt1) 
COMMON PZEROsHZERO 4 TZERO 
P=PZERO-XOUT(J) 

PRATIO=P/P7ERO 

CALL ENT (PoeTeY(1l) 92) 
TRATIO=T/TZERO 

CALL FORRHO (PeT»RHO) 

CALL FOR RHO (PZERO,TZEROsRHOZER) 
Z=P/(T#*RHO*1- 8644) 
RRATIO=RHO/RHOZER 
G=(P=$PSAVE)*(T+TSAVE)/(T-TSAVE)/(P+PSAVE) 
GAMMA=G/(G-le) 
GR=Z*GAMMA/(le—-(1le-Z2)*GAMMA ) 
VEL=0.6 
W=eZ2e"(HZERO=Y( 1) )*10336*98067 
IF(W) 69697 

VEL=SORTF(W) 
H=Y¥(1)%*e0242048+142,60 

PSAVE=P 

PSAVEST 

Peat d= 2424.3 

PRINTS sPRATIOsTRATIOsRRATIO2P 5 Te GAMMA 9GReoHoeZ s VEL 


RETURN 
PRINT1sPRATIOeTRATIOsRRATIOsP 5 Te GAMMA 9GRoeHoZ s VEL 
FORMAT (1290H1 P/PO T7110 RHO/RHOO p 
T GAMMA EFF GAMMA H 2 VEL 
9/36X% 84H ATM DEG: Kk 
CAL/GM CM/SECs7/73F l2e4eFl2e5sFilzele2l lZeS9Fl2e24 
oF l2e4sF 1201) 
FORMAT (3F12e4sFl2e3sFle2els2Fle2e4sFl2e2sFl2e4eFl12e1) 
RETURN 
END 
0932 
TOTAL O03 5* 
PA Ele 
SUBROUTINE FORRHO {P»sT RHO) 
R=1-¢8644 
VAR¥T/P=08255+162657*P)/(201* 1) 35 333 
RHO=1./V 
RETURN 
END 
0007 
TOTAL O007* 
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H6. 


P IN ATM T IN DEG K H IN ATM=CC/GM 
SUBROUTINE ENT (PoeTosHoL) 
DIMENSION HZERO( 130) eIHt130) 
HOFPF (PoT)=4e333*(82,6*#P+1-¢265*#P#P)/( eO1*T)*#*¥3- 333 
GO TO (1le2s3e4)eL 
4 J=C 
ti A at 
ZERO -INTBTUsEB-MOLE TH: IN VDEG oR 
READS »sHZERO( J) so TH J) 
3) FORMAT (2F10.0) 
aC = Pee ieG 
RZERCO(JI=HZEROI J) *2. 52081 
TPUTHE 69697 
NOW HZERO IS IN ATM=CC/GM AND T IS IN DEG K 
ZERO H AT O DEG GASEOUS 


6 JEND=J-1 
it RETURN 
2 00-339 J=leJEND 


LE Cl=TH O97) 32533429 
2) 6) CONTINUE 
GO TO 100 
32 TBeJ=] 
H=HZERO(I1B)+(HZERO( 1B+1)-HZERO(1B))/( THC IB+1)-TH( IB) )*(T<“THCIB)) 
1 -#HOFPF(P,T) 
RETURN 
2 DO 23 J=lesJENDs5 
HI=HZERO(J)=HOFPF (P,TH( J) ) 
LE Chl=HA) 23925024 
Zo CONTINUE 
J=VJEND 
24 JH=J=-5 
LP(SAi ss TO0s LO0Cy 29 
29 DO 27 K=JHeJ 
HA=HZERO(K)—-HOFPF(P »TH(K)) 
IF (HAH) 27928928 
27 HB=HA 
GO eu 
28 IB=K-]1 
T=TH(IB)+( THC IB4+1)—-THC IB) )* (H~-HB) /( HA~HB) 
RETURN 
Zo T=TH(J) 
RETURN 
100 PRINT1019PoToH 
101 FORMAT “742H TOR Hi 1S. OUT OF THE FRANGE “OF HZERO INPUT 
Lease 2k 02H Pel exe 2h lel 2% 62H Ha/Sh1 462) 
CAUL< EXIT 
END 
0047 


TOTAL 004 7* 
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ip Vitis 


MAIN TO PUNCH ENTHALPY TABLES 
DIMENSION HH 6) 
CALL ENT(P5ToHo4) 
T1=30.6 
DO 2 L=194 
T1=T1+606 
T2211 10% 
PUNCH7»9T2 
FORMAT (F17el) 
P=0-. 
DO 2 J=1%925 
P=P+.5 
T=Tl 
DO 6 K=1s6 
T=T+10¢ 
CALL ENT (Po9ToHH(K) 93) 
HH(K )=HH(K )*#20242048+142.60 
H TO BE PUNCHE D TO AGREE WITH DIN [ E GPC AIEZGM 
PUNCH3 »P sHH 
FORMAT (F6e1ls6Flle2) 
CALEY EX ET 
END 


HT. 


QO HAT GO SOLID 


TOTAL 
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H8. 
PART 3 


LABEL 
SUBROUTINE DIFFEQ (NoXeYsDY) 
DIMENSION Y(2)sDY(2) 
COMMON PZERO 
Sele 
P=X 
GO TO (192) oN 
P=PZ7ERO—xX 
C=—] » 
CALE ENT (Pele Y tl) +2] 
CALL FOR RHO (P.+eTsRHO) 
OY Cliy=S57RHO 
DY(2)=le 
RETURN 
FND 
OOT> 


TOTAL 0015* 


MAIN FOR SHOCK CALCULATIONS 
CALL ENT (PoT9Ho4) 
PRINT 
FORMAT (120H TOLL POl TO Pi tA 
1 P2 b PO2 TO2 PO2/P01 P1/P02 
2 9) 
CALL START(1) 
END 
0008 


TOTAL 0008* 


LABEL 
SUBROUTINE EQ1(V»REM) 
COMMON PZEROsTZEROsHZEROsTOLL 9 P29T29P19T1 9H1 
CALL FOR RHO (PleTisk1i) 
U1l=SQRTF(2¢*(HZERO-H1 ) ) 
U2=U1*RISV 
Pe=PitR i ful sUl=VA02e U2 
H2=HZERO-U2*U2* 65 
CALL ENT(P297T20H2Z92) 
CALL FOR RHO (P29T25R2) 
REM=R2-V 
RETURN 
FND 
O013 


TOTAL 0013* 


H9. 


PABEL 
SUBROUTINE ENT (PsTtsHol) 
PIN CAIM TINO DEG: Kk «oN ATM=e€C7Gm 
DIMENSION HZERO(130).»TH( 130) 
HOPPE (P67 34533934162 .68P +1 265*P FP) 7 (OT) #4396333 
GOV TO Ai s25354 304 
g=0) 
J=J+1 
HZEROMIN BIUsER MOLE ‘TH INJDEG R 
READS sHZERO(J)sTH(J) 
5 FORMAT (2F19.0) 
TH(J)=TH(J)/1¢8 
HZERO( J) =HZERO(J)*.52081 
IFC(TH(J)) 60657 
NOW HZERO IS IN ATM=CC/GM AND T [IS IN DEG K 
ZERO H AT O DEG GASEOUS 


“If 


A JFND2eJ=-]1 
1 RETURN 
3 DO 33 J=lsJEND 


Pe Chat) 3245335535 
33 CONTINUE 
GO 10-700 
a7 TR=J-1 
H=HZERO(1B)+(HZERO(IB4+1)-HZERO( IB) )/(TH( IB+1)—-TH( 1B) )*(T-TH( ITB) ) 
1 -HOFPF(P 4T) 
RETURN 
2 DO: 23° Ja leJIENDe 5 
HP=HZERO(JI—HORCr tes aio 
IF(HI-H) 23925924 
23 CONTINUE 
J=JFNN 
24 JH=J=5 
LPtghy) 1006 100629 
29 DO 27 K=JHeJ 
HA=HZERO(K)—-HOFPF(P,TH(K)) 
IF (HA-H) 27928928 
Tage HR=HA 
EO FO. 0 
22 gh | 
T=TH(IB)+(TH(IB+1)~TH( 1B) )*( H-HB) /( HA-HB) 
RETURN 
25 THTHt) 
RETURN 
100 PRINTIN1s»PoToH 
Or: FORMAT (742H T OR H [1S OUT OF THE RANGE OF HZERO INPUT 
Urey LAX 62H Pel2xXs2H Ts l2xXs2H Me7 3h i4s2)} 
CALLER ET 
END 
0048 


TOTAL 0048 * 





i 


SUBROUTINE PRINT (No gXOUTsYsDYo9J) 
DIMENSION Y¥(2)9DY(2)eXOUT(1) »R(3)sREM(3) 
COMMON PZEROsTZEROsHZERO,TOLL 9 P29T259P1l5T1lsH1 
GO TO (394) 5N 
PlHeP7ZERO=x OUT) 
CAEL ENT “(Ril bits yt Ll y<2) 
move yn) 
FANDUP ER ST. GUESS AlMRHOZ 
CALL FOR RHOUP Le lis2(1)) 
CALL FOR RHO(PZERO¥TZEROSRZERO) 
R(1)=1-01*R(1) 
R(2) = RCT) 
DR=el*(RZERO-R(1)) 
CALL EQ] (R(1)sREM(1)) 
R(2)-="RUZ) + DR 
CALL EQ] (R(2)sREM(2)) 
TEC REMCITY*REMUZ) ) 106 10511 
RiS)J=CR( 1 AR 2) te 5 


CAELVEG] (R(3)-6REM(3) ) 

L=1 

Pe c6R EM 2 REM ol 261 56 14 

L=2 

RE yaR tS) 

PE CCR( 2) Ri) eee 7 R= TOL es a6 10 
CONTINUE 


CALE STAR I C2) 

TF CY Cl) HHZERO}) 33656 

XB=Y(Z2) 

HB=Y(1) 

RETURN 

POZ=¥i 2) (Yt 2) =X BAe yl Py KHZERO) IY Sie) 
CALL ENT (PO2sTO2ZsHZEROs2) 

POUT=PO02/PZFRO 

POUTA=P IVP Oe 

PRINT? sTOLEsPZEROsIZEROsP les Tllseee2sT2ePO2stC2sFP OUT sPOUTA 
FORMAT (F106694(F10.39F1001)92F1065) 

CALE “SIAR 1 1) 

END 


TOTAL 


H 10. 


0038 


0038* 





1 


SUBROUTINE PRINT (N pXOUTsYsDYoJ) 
DIMENSION Y(2)+s+0Y(2)sXOUT(1) ¢R(3)eREME3) 
COMMON PZERO,TZEROsHZEROs TOLL sP2eT25P1l9T1l Hl 
GO TO (394) 5N 
Pl=PZ2FRO-XOUT(J) 
CALL ENT (PloT1lsY(1)92) 
H1=Y(1) 
FIND FIRST GUESS AT RHO2 
CALL FVPOR RHOCP is 713R 0 15) 
CALL FOR RHO(PZEROsTZEROsRZERO) 
R(1)21-01*R(1) 
Rete) = tC 35) 
DR=.1*(RZERO=R(1i)) 
CALL EOW. CRC s+ REMC) ) 
Ri2) =] R(2) +DR 
CALL EO] AR(ZISREMCZ)) 
[FP CREM( 1) *REMIZ2)) 1] 0610511 
RUSyY=( RC 1)ER U2) ) 25d 
CALL EQ1 (R(3)sREM(3) ) 
L=1 
IF (REM(3)*REM(1)) 12913513 
ie 
R(L)=R(3) 
IFCCR(2)2R(01))*2e¢/R(1)-TOLL) 14914510 
CONT IT NUE 
CAUE START (2) 
IF (YC1)“HZERO) 52655 
XB=Y(Z) 
HB=Y(1) 
RETURN 
POZASY C2) = ty (2 y=x 8) CY (ly -HZEROY7 CY Gl) -i6) 
CALL ENT (PO2sTOZ2sHZEROe2) 
POUT=P02/PZFRO 
POUTA=P1/P02 
PRINT7s TOLL sPZEROsTZEROsPloeTlsP2sT29P029T02 POUT sPOUTA 
FORMAT (F106.654(F 10, 36F 106 l)s2F 1065) 
CALL START 14 
END 


TOTAL 


H10. 


0038 


0038* 





100 


200 


SUBROUTINE START(N) 
DIMENSION MARK(4) »XOUT (200) »Y(2)e9TOL(2) sYMIN(2) sDY(2) 
COMMON PZEROsTZEROeHZEROs TOLL »P25T2 


H= 095 
DO a THA 9:2 
MARK (J)=1 


YMIN(J)=20¢ 

GQ TOt(100%9200) oN 

READ4G, PZERO»sTZERO»eP1»TOLL 

SET CONSTANTS FOR RUNGE INTEGRATION DOWN THE [SENTROPE 
X=06 

TOW Ty TOL 

CALE."ENT) CPZ EROS IZEROeY Cli s3) 
HZERO=Y (1) 

MARK (4) =0 

FORMAT (4F12.-0) 

XOUT(1)=PZERO—={P1 

CALL RUNGE (19Xe¥e2TOLsYMINoHs XOUT »MARK) 
CAI ExT 

X=P2 

CALL ENT €P2eT2sY¥tlis3) 
MARK(2)=(PZERO=jP2 )*200e+206¢ 


TOL(2)=le 
Y(2)=xX 
XOUT(1)=P2 


MAR=MARK (2) 

DO 7 K=2sMAR 

XOUT(K)=XOUT(K-1)+e905 

CALL RUNGE (29XsYsTOLsYMINesH»sXOUT »MARK ) 
CAL Est 

END 


TOTAL 


SUBROUTINE FORRHO (P»T»RHO) 

R=] 2.8644 

VeR PIP = (626 tle 2 ODA Hy, Ue Cll) e% 36533 
RHO=1e/V 

RFTURN 

FND 


TOTAL 


H1i. 
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PART 4 


LABEL 

PROGRAM FOR RATE CONTOURS 

READ SECTION 

DIMENSION R(5)9 PP(5)s FP(3)s P(3) 
READ1s RHOAsRHOB»SURFA»sSURFB»WMOL 
FORMAT (5F12-0) 

READ1» TSTARTsTENDsDELTAT 

READ4O »R 

FORMAT (5E12e2) 

PRINT22s RHOAsRHOB »SURFAsSURFBsWMOL 5R 
FORMAT (1L4HIRATE CONTOURS »// 


6H RHOL=9F8e3s 2H —s F9e6s10H (T DEG C)o/ 
8H SIGMA =9F8e392H —sF9e6910H (T DEG C)o/ 
THEWMOL =6F Ge ls 77 
24H TABLE CONTAINS P IN ATMo// 
19X96H RATE =0E612¢394E18632s10H P INF o/ 
SA. 1 DEG. K 7 

5H =9/ ) 


M=(TEND=TSTART)/DELTAT 
T=TSTART=DELTAT 

DO: 20° [=lem 

T=T+DELTAT 

COMPUTATION 

RHOL =RHOA-RHOB* (T-273e16) 

SIGMA= SURFA-SURFB*(1T-273216) 

CALL PINFIN (TsPINF) 

IF (PINF) 205920574 

C= eOlLS25ES4 14061671 1.636% |) FP INE 
C1=LOGF (C1) 

C L=LOGF ( SQRTF(CWMOL¥2e¢*SIGMA/ (3461416%*64025E23) )/RHOL )+2e*Cl 
C2=2e*S1TGMA*®WMOL/(322062*1201325E6*RHOL*T ) 
C2=C2¥C2*4,18B8*¥SIGMA#120E16/(1e38*T) 
NO 19 J=195 

IF(R(J)) 19919921 

A=LOGF(R(J))=-C1 

B=—26 

C=C2 

CHOP FOR P 

P(3)206 

FP(1)2C 

P(3)2tP(3)+36 

IF(P(3)-554) 66966319 

IFC FP(1)*(P(3)*#8(3)* (A4B*¥P (3) )4+C))2910011 
P(1)=P(3)-3e 

P(2)=P(3) 

P(3)2=(P(1)4+P(2))%*e5 

IF (P(2)-P(1)—-2e005*P(2)) 1091099 

L=1 

PPI3)=P 03) #P (3) 8 (AFB FP (3) 4G 

IF (FP(1)*FP(3)) 491095 

L=2 

P(L)=P(3) 

GO TO 6 

IF (P(3)=-70¢4) 18518919 
PP(J)=EXPF(P(3))*#PINE 





- i 
Die 
| 


~ 


“ 


CONTINUE 

PRINT =32% 
FORMAT ( 

GO.710 -2101 
END 


Ty PP, PINF 
F10.299X95E18e3eF1l163 


) 


His. 


FORA 


0060 


0060# 


fi 





100 


— 


15 


24 


81 


25 


PART. 5 Hoke 
MAIN FOR CONDENSATION 


DIMENSION Y(10)+TOL(10)sYMIN( 10) »sXOUT( 100) sMARK(5) 
COMMON PZERO»s TZEROe+WMOL » GAMMA>s RHOL » 


L1HFG»s SURFA, SURFBse XIe ALPHA» AMACHsy Xo AREA, 
2ZDADXs Ps Ts RSTARs TDs AMASS»s AVGR» ASTARs 
3 SSAT,» INTERP eZ »sYDROPsYRATE » GAMA »GAMB oY sNPART 


»sRTIMES 
INTERP=1 
CALL AREAC(XsAREAsDADX) 
READ 1» PZERO ATM 
READ 1» TZERO DEGe, GC 
READ 1» WMOL 
READ l» Z 
READ 1+ GAMA 
READ 1s GAMB 
READ 1, RHOL GM/CC 
READ 1» HFG BTU/LBM 
READ 1» SURFA DYNES/CM 
READ 1s SURFB DYN/CM=C 
READ 1l» 
READ 1l»e ALPHA 
READ 1s RTIMES 
READ le TOL(4) 
READ 1% TOE 
READ 1+ AMACH 
READ 3+ CNTAM PTS/CM3 
READ 35 RZERO CM 
FORMAT (€F 2060) 
FORMAT (Elle3)} 
PRINT2s PZERO»s TZEROs WMOL >» 2+GAMA »sGAMB»sRHOL »HF Gs SURFA 


l1»sSURFBs XI» ALPHAs RTIMES» TOL(4)5 TOLL» AMACHs CNTAMs RZERO 


FORMAT (25H1 CONDENSATION IN NOZZLES» ///9H PZERO = 94 
F1l0e3%9/9H TZERO = 5» FlQNe2s/7HWMOL = 9» Fl2e29/5H Z = 9F1l2e3s/7H 


2 GAMA =9Fl1l2e29/9H GAMB = »9F1l2e¢6e/8H RHOL = oFl2e39/7H HFG = -» 

SP l2e2z97 7H! SUREA =) 9.) F1Os24e/9H SURFS = & Elless/6H XI = » Fl6e5¢7 
48H ALPHA =9F14e59/11H RATE*100** oF8e29/9H TOL(4) =sF18e8s/7H TOLL 
5=09F18¢8s/9H AMACH] =9F10e2e/8H CNTAM =9E1l2039/8H RZERO =sEl2e3s//) 


PZERO=2115-PZERO 
TZERO2#1¢8* (TZERO+273-16) 
INTERP=2 

XOUT(1l)=le 

XOUT(2)=AMACH *AMACH 
XOUT(3)=2¢*XOUT(2) 


H=0e0] 
MARK (1)=1 
MARK (2) =39 
MARK (4) =0 


DO 24 J[=599 
VOCs =7 CL 

NO 81 T=193 

TOL (1)=0¢1 

DO 25 [=199 

YC = 05 
YMIN(1)=9200001 





WMOL=WMOL /2 
¥(1)22501327*6204*CNTAM* 1544.*TZERO/( WMOL*PZERO) 
Y(2)=Y(1)*RZERO 
Y(3)=Y(2)*RZERO/26 
Y(5)=le 
Y(6)=1e 
P ZERO=PZERO/2115-¢ 
TZERO=TZERO/ 108-2736 16 
NPART=1 
X=O0-6 
N=9 
CALL RUNGE (Ns X* Ys TOL» YMIN»s He XOUT» MARK) 
READ 779 L 
Te FORMAT (12) 
Pe Wie ON 1015 100 
101 CALL EXIT 
END 


TOTAL 


* yea Ib 
C SUBROUTINF TO COMPUTE NUCLEATION RATE 
SUBROUTINE RATEC(PRESSeTEMPs SIGMAs SSATs RATE) 
DIMENSION Yt10) 
COMMON PZERO»s TZERO»WMOL»s GAMMA»s RHOL » 
lHFGs SURFA,», SURFBs XI» ALPHA» AMACHs, Xs AREA>s 
Z2DADXs Pe Teo RSTARs TDs AMASSs AVGRe ASTAR»s SSAT 
3 »sINTERP 92 sYDROP s YRATE » GAMA 9 GAMB 9s YoNPART »RTIMES 
A=RSTAR 
B=10.**16/(1.380*TEMP) 
C=-4e1)887*STGMA#¥A#A#B 
3 D=PRESS#1.01325E6*B 
14 E=WMOL/ (RHOL*#6e925E23) 
FEZ2e*#S IT GMA/(3e¢1416*WMOL ) *62025E23 
AA=2e*LOGF (DI LOGFILE) +eo*LOGFIEEI+ © 
AA=AA+2¢303#1 66 
TFCAA) 16916417 
16 RATE=0,.0 
GO TO 12 
ye CONTINUE 
IF (AA=8520) 18918919 
19 AA=85,0 
18 RATE=EXPF(AA) 
l2Z RETURN 
END 


TOTAL 


H15. 


0072 


O072* 


O0i25 


O10i2 3% 





H16. 


LARFL 

SUBROUTINE TO COMPUTE DERIVATIVES 

SUBROUTINE? DPEFEEG @UNs xs. Ys DY} 

DIMENSION “¥ CTO) sDY 010) 

COMMON PZERO»s TZEROysWMOL » GAMMAs RHOL » 
1HFG»s SURFA,» SURFBs XI» ALPHAs AMACHy Xs AREA, 
2ZDADX» Ps Ts RSTARs TDs AMASS»s AVGR» ASTARs 
3 SSAT, INTERPs2 »YDROPsYRATE »GAMA 9GAMB sY eNPART 

ao TEMP=Y(6)*(TZERO+273416) 
=,OR847B/WMCL 

GAMMA =GAMA—-GAMB# (TEMP-273216) 

FACZHLle/(let(GAMMA-1le)*(1le-Z)) 

GAMMA =GAMMA*F AC 

GO TO (192) sNPART 

DO-3> 1=1s59 

3 DY Cia =0% 

DY(5y2-Y(5)/ 020 /GAMMA+X* (GAMMA-1.¢)/GAMMA ) 

DYUCBJ=RSTEMPEDY (577 YU5) 

DY(6)=Y (6) *(GAMMA-1.)*DY(5)/(GAMMA#Y(5) ) 

RETURN 

2 CALL PINFIN (TEMPosPINF) 

PRESS=Y{(5)*PZERO 

SSAT=PRESS/PINF 

Ler ASS A aie) 205° SOs 231 

30 RATE=0 

RSTAR=0 

DY(9)=06 

DROPG=0 

C12 UO mae 

a1 SIGMA=SURFA=-SURFB* ( TEMP-273-216) 

RSTAR=2e*SITGMA*WMOL /(824062*1.201325E6*RHOL*TEMP*LOGF(SSAT) ) 
IF (Y¥(1)=25.1327) 655965966 


~ 


65 AVGR=1.23*RSTAR 
SO 10°67 
66 AVGR= SOR IT U2e* YC S)7y¥ (1) ) 
67 CALL RATEC(PRESS» TEMP» SIGMA»s SSATs RATE) 


CALL DROPGC{PRESSs TEMPs SSAT»s AVGR» DROPG) 
DROPG=DROPG/(Y(7)*SQRTF (GAMMA#8e314E7/WMOL# TEMP ) ) 
VEL=SQRTF (—-2e#Y(8)*980-7E6) 
DYL9Y=E1le /VEL 
a2 CALL AREAC(X» AREAs DADX) 
CONLAM= (GAMMA=120)/GAMMA#HFG#WMOL/(1544,#1e¢8* TEMP)*778,4 
DY(1)2=25e1327% RATE#AREA 
DY(2)=Y (1) *DROPG+DY(1)*RSTAR 
DY(3)=Y(2)#DROPG+12.56637* RATE*AREA#RSTAR#®RSTAR 
YDROP=RHOL /AMASS#Y(3)*DROPG 
YRATE=RHOL /AMASS¥®#4e¢1888#RATEXAREARRSTAR## 3 
DY €4 )=YDROP+YRATE 
DO 444 J=1 94 
444 DYt1-4 20 
DY(5)=Y(5)*¥( CC (CONLAM-1e/(1le-Y(4)) )*DY(4)—DADX)/¢t(1le (le-Y(4))* 
1 ( (GAMMA=1¢)/GAMMA+1¢/(GAMMA#Y(7)#Y(7)) ))) 
DY(6)=Y(6)# ( CONLAM#DY(4)+(GAMMA—-1¢)/GAMMA#(1e~Y(4)) *DY(5)/Y(5)) 
OY Cy =¥ CPA HS Hy (4) VED Y ES) 7 UGAMMASY ET yey (7 ey 
] e27Y C6 )*DY16) ) 





(ns Ten SD Fe 2. 


al 


10 


R1T.! 


DY(8)=R*TEMP#DY(5)/Y(5) 
RETURN 
END 


PINFIN FOR CO2 YW TO 240 DEG K 

BELOW SS Pl NPaCeE RPE Centos) Wie -USeDUANDUP TINE Wiles - USEO Al “Gee 
ABOVE 93 AN EGUATIUN OF THe FUR PEFCHEXPF(A/T) wAS Used TU [NWT eERKPULATE 
BETWEEN DATA POINTS. DATA TAKEN FROM HAUDDUUR UF RPHYSTCdS AND 
CHEMISTRY» INTERNATTONAL ZRITICAL TABLES» AND MATHESUN GAS DATA 


SUBROUTINE PINFIN (ToPINF) 
ie (T-93e¢)lols2 
A=6.VU33 
B=-1263. | 
C16 in Oe Ae. 
IP (114606) 39394 
A=72492 
B=-1430-6 
GO O10 
IF (T1600) 55596 | 
A=7.0851 
Se-13736 
GO. 20x20 
IF (T=195¢e) Vo798 
R=6,.6 555 
B=-1333-. 
GOr 10: yo 
TP A 1=216060) Sse 11 
Az=7.028 
D=-1367-6 
GO .10 20 
Az4,72?7 
E z—-Y70el 
FINF=1lOe**(A+B/T ) 
RETURN 
END 
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26 


27 
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67 


68 
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LABEL 
SUBROUTINE TO PRINT 


SUBROUTINE PRINT(Ns XOUTs YOUTs DY» J) 


DIMENSTON =xXOUTU])s YOUT C7} 9Y(10) 
COMMON PZEROs TZERO»sWMOLs GAMMAs 


LHFGs SURFA» SURFBs XI» ALPHA, 
ZDADXs Ps Ts RSTARs TOs AMASS >» 


Pet d=2) ec2le22e 32 


AMACH, 


sDY(1) 
RHOL » 


Xe AREA, 
AVGR»s ASTAR» 

3 SSATs INTERPsZ sYOROPsYRATE »GAMA»GAMB sY sNPART 
GAMMA=GAMA—-GAMB*#YOUT (6)* (TZER0O4+273016) 
GAMMA=GAMMA/( 1le+(GAMMA=1le)*( 1le-Z)) 


HG; 


AMASS=ASTAR#¥PZERO*YOUT (5) *2115¢*SQRTF ( GAMMA¥WMOL /(3202*15446 


le SAT ZERO+27 Je 16) xOUT(S) 153 )) 


AMASS=3262* 45406 * AMASS/ (30 048% 30448) 


PRINT1O 
PSTAR=YOUT (5) 
TSTAR=YOUT(6) 
GSTAR=GAMMA 
RETURN 
NPART=2 

DO 26 IT=199 
YU) -yYouT.b) 
DO 27 I=193 
Y(OT)=YCT)*AMASS 
Y(7)=AMACH 


A=PSTAR/(Y(5)*AMACH )*SQRTFCY(6)*GSTAR/(TSTAR*®GAMMA) ) 


INTERP=3 

ARFA=A*ASTAR 

CALL AREAC (X»sAREAsDADX ) 
INTERP=2 

XOUT(3) =X 

DO 25 12363 
MOUT(3E+1)=XOUT( 1 )+e2 


CALL AREAC(XOUT(J)»s AREAs DADX) 


AREA=AREA/ASTAR 
TEMP=YOUT(6)*(TZERO+273016) 
CALL PINFIN ( TEMP» PINF) 
PRESS=7OUT C5 1 *P ZERO 
SSAT=PRESS/PINF 


S IGMA=SURF A=SURF B*# (TEMP-273016) 


DEC SSAT=16) 69969967 


CALL) RATEC(PRESSs TEMPs SIGMAs  SSATs 


PP CYOUT (IT )=25e1327) 695694568 
AVGR= SQRTF(2e¢*YOUT(3)/YOUT(1) 
CALL DROPGC(PRESSs TEMPs SSATs 
DIDETO=Y OUTS) UTZEROF27 3616) 
PPSF=2115e¢*PZERO 
TDOEGR=12¢8*(TZERO+273416) 


TOTALN=YOUT(1)*WMOL *PPSF/(2541327*62 e4* AMASS*®1544e*TDEGR) 
PRINT11XOUT(J) sAREAsYOUT(S) s YOUT(6) sSSATs YOUT(7) sYOUT(4) » 


) 
AVGR >» 


RAT 


1 YDROP sYRATE sDTDsAVGReRSTARsRATE s TOTALN 


RETURN 

FORMAT (120H x A/A* 
1Y GROWTH» NEW PART TO-T 
2 //) 


P/PO 
AVGR 


erg) 


E} 


DROPG) 


S SAT 
RSTAR 


MACH 
NRATE 


MOIST BY GROWTH 
TOTALN 





el FORMAT 
END 


Le) 


10 


i 


60 


H19. 
2 


(Fl eZeFBe3s2FleGsFBe SF fe SeFB ede ZELOeZsFle2v4FE10—3) 
0057 
TOTAL 0057%* 
LABEL 
SUBROUTINE DROPGC{(PRESS»TEMPs SSATs AVGR»s DROPG) 


SUBROUTINE TO COMPUTE DRUP GROWTH 
COMMON PZERO, TZEROsWMOLV»sGAMMAV»RHOL » 


IHF Gs SURFAs SURFBs XI» ALPHAs AMACHs X»5 AREA, 
2DADX*+ Ps Ts RSTARs TDs AMASS»s AVGRs ASTARs» 
2SS5A0 eT hTERP 


J=1 

P=PRESS*1.01325E6 

f=TeMP 

R=AVGR 

L=0 

D¢«1)=0 

D(2)=0 

TD=T-le 

PDS TD ral. 

CALL PINFIN(TD»s PINF) 

PD=tl-eOlS2566*7 INFFEXPECROTAR* It TOFRy) SLOG CSSAT 
SASCON=E2314E7/WMOL 

CONL AM= ( GAMMA~1.) /GAMMA#HFG*778./(1544.45*1 ed*7 SMP) #WMOL 
BETA=P/SQRTF (662832*GASCON*T ) 

BETAO 2 -PD/SGRIE (6.28227 GASCONF DO) 
DROPG=XI/RHOL*(BETA=BETAD) 


AK = le - BETAD#TD/(BETAT} 
AAA= e5* (CONLAM=1 « ) ®GAMIMAS( GAMMA-1.)*(1e-BETAD/BETA) 
A= AA+tAAA = (1e/XI —- le) #¥ALPHA*(TD/T=-1.) 
GO TO (839) s5J 
wae 
IFCA) 10381099 
DROPG=O-64 


PRINT11» PDsPsTsAsAAsAAA »XeV¥ (4) 
FORMAT(SE12.3 ) 

et RN 

T=1 

IF (A) 29293 

T=2 

D(I)=TD 

DS] Drie 2) 

IF (D3) 459194 
LB=.2*1D0CE)+pt2 7) 

L=L+l 

LP tli=20)) 3455.60 
IP UDRORG) WOs5 36 

RETURN 

END 





re 


ONOOD * 


21-0 


220 


230 
250 


1030 


1040 


H 20. 
LABEL 


PIRSTSORDER OLEPs-eG. ROU ING =-ADJUsS sae teEPo Sore 


SUBROUTINE RUNGE(Ns Xs Ys TOLs YMINy Hs XOUTs MARK) 

DIMENS PON: “YC 1) 5 YANO eo ROE CT Ys SUBL SO lis XOUT U1 jes rian th) 
DIMENSION” (DY CSOs 1A (30 ls iF RESO)6 Bt S50) se (FC (50)'s. YREEP (50) 
KBTWN = 1 


NCOUNT = 15 

J = MARK(1) 

MAX = MARK (2) 

L = MARK (4) 

[key Abies) 2g oc Oss e 20 
RE Sie. v1 


OOa 2505 desis iN 

SUB CIS) ceo WO 3:20 
PEM AX I, 205 305. 30 
RETURN 

A = XOUT(uU) = X 

B = ABSF(1.E-7 * X) 

he KGa 3) “Os 32455 33) 
IF (A - 3) Ds 50s 60 


OP ia ee ar 

GO-TO <0 

CALL PRINT(Ns XOUTs Yo DYo J) 
a ane ae oe t 

GO6TerieG 

hE-1GA. = pte S 10s," (Os Go 

H = A 

Lea ae] 

oO. [Oo qe 

IF (A = 3¢%*H) 90+ 1000» 10090 
H = eD*A 


DO RUNGE=RUTTE=MERSON INTEGEATIGON 


XA XP OAS 3% 
XG = aa 0 5%*H 
COULADU PEREGO UNS BX 542) 6 1D} 


Ht 


X = X + H 

DOO 30 I = lo N 

VREEP CLI = Fit) 

FAC]. = A*DY (1) 

YACIG 4] YCI)<+°FA(TI73. 


CALL’ DIFPEQUNs -XAs Yas DT) 

DO 1040 | v=o 16- N 

YAGI) p= VOT bt PACT 60> * Web (13766 
CALL DIFFEQ(Ns XAs YAs DY) 

DO-.1C 59 I = ls N 

BS) Gan Oo 


%) 





1030 


1060 


1090 


1190 
ep 
Ea 8; 
135 


le Rete 


1140 


1150 


1160 


1i6s 


LO 


Lii2 


1174 


iio 


1178 


Hel. 


YA eee Woy ites leo A Gl 4 6 S37 ORF EIT 
CALLS IDIFFECUNs  XBis. YAs. DT} 

DO 1060 I = ls N 

Pet) = (repre) 

YO a= erate oer A Ieee ctr et lit 7 tec] 
CALL DIFFEQ(Ns Xs YAs DY) 

DOs 130 I = l»s N 

Y(I) = YCI) + FACT) 7/606 + e666666667*F C(I) + H*¥DY(1)/66 
Cee a On 

IF (ABSF(U) - YMINCI)) 12305 1090s 1090 
KLOW = 2 

Eno 22 AGSr COR —4 YA 

LPM ee ceva BoE CIOL CL ical). ) LET Oerl LOOs L100 
KBIG = z 

GO PCO 4.36 

le a Je = ASSetSUBA 1) eu). ) PS O607 ZO se bic 
KBTWN = 2 

CONTINUE 

GO TO (1004, 1135) + KLOW 

GO TO (1180s 1140)» KBIG 

NCOUNT = NCOUNT - 1 

be MI COUN TY 2 S0s 0) 50 ssl 70 

PRINT 1160s Xs H 

PRINT -P165s “(ler YU). DYE); I = ls N) 


RETURN 

FORMAT 156HS4S1EP SIZE HALVED 15 TIMES. CONSECUTIVELY “SINCESEAST {PRI 
INT /29H PROGRAM TERMINATED AT X = s El6s8s- SH H = » El6eds 
Z/7 3M ] 9 13Xs SHY(IT)s 16Xs SHDY( I) »//) 

FORMAT (13% 7X8 2(E1668s 4X)) 

KBIG = l 

Les OCR = By Tacs el 2 sly 

X = X = H 

H = «D%#*H 

DO 1174 I = ls N 

Y(1) = YREEP( 1) 

KBTWN = 1 

KLOW = 1 

Go 70. 2000 

M= 15 = NCOUNT 

PRINT 1178s Ms Xs H 

PRINT T1655 (is YC1]l)s: DYlI)s 1-= Leon) 

RETURN 

FORMAT (41H4STEP SIZE BECAME TOU SMALL FUR COMPUTEKs/20H IT HAS GE 
1EN HALVED » [2s 21H TIMES CONSECUTIVELY e/29H PROGRAM TERMINATED AT 
2 xX = 7 E1660 46 ~Ohis H = » El665s//53H Is 13X99 GSHY( I)» 16X»s 
SSADY RT IS 7) 

NCOUNT = 15 

SO TO (1190s 1200).s KBTWN 

H = 2e%#H 

KBTWN = 1 

KLOW = 1] 

CHEEK FOR INTERMEDIATE - PRINT OUT 
GO. TO C106 =10)5" (FIIEST 





1050 


1060 


1090 


1190 


AO 
wae) 
L130 


35 


1140 


DES 0 


1160 


Nid Efe) 


1170 


Lie 


1174 


slap aye) 


ee ys) 


H2l. 


WA = ale ee el ZORA + God or ro) 

CALL DIFFEQIN»s XBs YAs DY} 

DO O60 I = lo N 

FC(1l) = HeDyY( 1) 

VAG) = YC oe Str Al 1) = ties ktr oe ( [to ether) 

GALE OLFFEOGNs, XY Aw Dy) 

DO i136 I = ls N 

Y(I) = Yl) + FAI) 766 + «66666666 7*FC( I) + H*¥DY(1)/6- 
J = Y(I) 

IF (ABSF(U) = YMIN(I)) 1130s 1090+ 1090 

KLOW = 2 

Eras Z*ACSEUU.—- YAL1) ~) 

Le St Be AB SEULOLGiIa os 2 LerOs 211C00%; -E1Go 

Re Gers =2 

GOGO) 1.130 

boo Base Ae Se ESUB CT ye) 3 Pi3s0's 11209) -Li20 

KBTWN = 2 

CONTINUE 

GO-1O, C1004. 41 135% “KLOW 

GO TO tll80s5 1140)s KBIG 

NCOUNT = NCOUNT - 1 

Feo COUNT \) elses] SOs 170 

PRINT 11609 Xo H 

Pin! W636 0 (cleat Cy eo yo) I = ls» N) 

RETURN 

FORMA <hSoHeSTEr SIZE HALVES 25 -TIMeEs CONSECUTIVELY “SINGE CLAS T oP 
INT /29H PROGRAM TERMINATED AT X = » El6e8s 8Hs Ho= =" El66ts 
27 FSH [Ie er Taxes SAYUT hs Tox SHDY (16773 

FORMAT (1359 7Xo 2(E16e8s 4X)) 

KBI1G. = 2 

Te (GH =) Ply. Wit 2e Lie 

X = X -H 

A = goth 

vO L1t4 I = ls N 

Yt) ic] “YREEP CT) 

KBTWN = 1 

KLOW = 1 

GO; TO 1000 

M = 15 = NCOUMNT 

PRINT 1178s Ms Xs H 

PR OUNY WO oe ot len eC ys: DY lies a ee les 

RETURN 

FORMAT (4l1H4STCP SIZE BECAME TOO SMALL FUR COMPUTEK./20H IT HAS BE 
LEN > HAEVED. +. 112s 21H TIMES CONSECUTIVELY s/29H "PROGRAM. TERMINATED =AT 
2 X = » El6.8s 8Hs H = 9 El6¢59//3H Ie 13X9 GSHY( I) s 16X»s 
35HDY(1)9//) 

NCOUNT = 15 

GSO TO (1190+ 1200)s KBTWN 

HH = 2e%*H 

KBTWN = 1 

KLOW = 1 

CHECK FOR INTERMEDIATE PRINT OUT 
GO Os (10s"_ 10) 5. (ETE St 





110 
120 


130 


140 
12.0 


PE tres) 120% 120». 
PRINT 130, Xs H 
PRINT 1405 (Is YC 
L = NUM 
GO" 10-10 


FORMAT (5H X = 4 
L1I6Xs SHC Yt) 7) 
FORMAT (55Xs 13>» 


L=Ll- 
GO TO 10 
END 


L 


150 


I)» 


El6. 


7X9 


DY Chis 


= l» 


Bs 4Xs 4GHH = 


2(E1668s 


4X) ) 


% 


N) 


El6e3>» 


11X» 


LHI» 


H 22. 


13X» 


GHY( I)» 





{\ 


at a ea, ae et as a ON Tc | 


10 


TZ 


16 


eg 


CY 


COA ACE. TA 


Vee) 


LABEL 
INTERPOLATION ROUTINE FOR COMPUTING AREA 

DAD LS CON TENUOUS 

SUBROUTINE AREAC( XX »AREA sDADX) 

COMMON PZERO»s TZEROs WMOL, GAMMA, RHOL» 

LHFGs SURFA»s SURFBs XI» ALPHAs AMACHsXX» AREA» 
ZDADX oP oe TesRSTARsTD»sAMASS »AVGRIOASTAR 9 

S3SSATs JATERP 

DIMENSION DDOUS0) ext 5D) A500 )sD150)sDA150) 


AREAC PART 1 


INTERP=1, READ AREA X AND ASTAR 

SEM ixsDACX. TO CORRESPOND )-T1OSAS TAR 
ASTAR MUST BE IN AREA DATA LEST 

GO TO (1s29394)5 INTERP 

1=0 

T=]+l 

READ1O, X(I)sA(1) 

FORMAT (ZF 15.0) 

Vea Tice Ves 2 1 

READ10O,s XX»sASTAR 

IEND=I-2 

PEND CS=—GAST VALUE CF T (USEDeEMINUS ST 

AREA=ASTAR 

DADX=0-6 

PING, ZERO 

I=0 

BES oat 

bet x Gl y= xx 132914914 

IZERO=] 

D(IZERO)=0-6 

DO 15 “T= I1ZEROs TEND 

D(T+1l1) = 2et*(ACT+H1)I—-ACT)I/(XC1L41)—-X01))-DC1) 

DO PG: (Fs=25 ZERO 

J=1ZERO=—1 +1 

Didyezet (ACOt At) 7 xt Xabi +1) 

DO. P= e2s LEND 


(0) Gia ea ae a sc ne ar aa?) Ae oa a a 
PEUG Sa eNO IBD 
RETURN 


AREAC PART 2 


GIVEN xX OUTPUT IS AREA AND DADX 
DE OXA= KCL 220 92 3422 

I=I-l 

lige wens fhe > Ce a ae ne 20%25%5.24 

IB=] 

LA=I+1] 

GO TO 24 


H 23. 





Hek, 


22 T=[+1l 
bee ewe 60) 254 23522 
25 [A=] 
IB=I-1 
é FIND AREA AND DADX 


24 AREA=ACIB) 4 (XX=X(16))* (DO( 14) 
ie OX ACB oo a (yD 
302 DADX=(9D(1A)*(XX-X(1I58))+0(16))/AREA 


RETURN 
c 
© 
C AREAC PART 3 
C 
c 
G GIVEN AREA FIND X AND DADLX 
C le MACH IS LESS THAN 1 LOOK BEFORE THROAT 
C IF MACH IS GREATER THAN 1 LOOK AFTER THROAT 
3 I=IZERO 
le CAMACHH=l0 3 20531551 
a0 T=]-1 
[TP tAREA=A17) 1) 32432950 
a2 IB=I 
TA=[+1 
GO TO 34 
cal T=I+1] 
IF (CAREA-A(T)) Sy ois | 
BMS I[A=I 
IB=I-1 
C FIND X AND DADX 
34 LP DUTANH=D.1B)) (365357525 
36 SIGN=1l. 
301 KX] (=DCIBJ+SIGN® SQGRIFU(DC( IE) #0 1B) =20*DD( 1A) *(A11B) =AREA))) 


be DOCVAY 2 -X4e084 

IF (XX—-X( 1B) ) 78939939 
38 SIGN=-le 

GCG" 10: 3401 
39 LE UAK=XC TAY e025 302538 
37 AXX=CAREA-=ACIB))/7O Ol). 4° A018) 

OADX=D(18)/AREA 
4 RETURN 

END 





a 


SUBROUTINE AREAC (XsAREA,DADX) 


Hid. 


COMMON PZEROsTZEROsWMOL »GAMMA »RHOL sHF Gs SURFAsSURFBsXI sALsAMeXXeAREA 


LAsDADX ePesTsRSTARe TDs AMASS es AVGReASTARSSSATsINTERP 
GO, TOstls2s3) s+ INTERP 
CONTINUE 
ASTAR=le 
X=e55 
RETURN 
NOZZLE TWO B 
IF (X=1e¢45)10011l911 
AREA=1e¢+e00889* (X—-255) 
DADX=e-00889/AREA 
RETURN 
IF (X=-4285)1291313 
AREA=12-008+.0218* (X=-1.45) 
DADX=.e0218/AREA 
RETURN 
AREA=12082+2e0133*(X—-4e85) 
DADX=.20133/AREA 
RETURN 
IF (AREA=-12008) 3023031 
X=(AREAwm1Le)*112¢5+055 
DADX=,.00889/AREA 
RETURN 
X= (AREA=12008)/20218+1245 
DADX=.0218/AREA 
RETURN 
END 


TOTAL 


0028 


0028* 





Tal. 
APPENDIX I 
ESTIMATION OF EXPERIMENTAL ERROR 
TEMPERATURE MEASUREMENTS 
Estimated error (difference between measured and actual temperature) 
on temperature readings except where otherwise noted is + Cereaet 1.5°F, 
This was determined after the following meaSurements and observations. 


1) Temperature of various portions of the ice bath of the reference 
junction were found to vary by as much as 1°F from the freezing point, 


32.2°F, despite its containment in an insulated thermos. 


2) A calibration of the potentiometer (Wheatstone bridge type) used for 
measurements, compared exactly with those readings of the standard copper- 
constantan conversion tables at each end (32°F and 212°F) of the cali- 


bration and showed a maximum variation of + 1.0°F between. 


3) The time constant of the thermocouple was measured with the aid of a 
Sanborn Recorder capable of recording a 50 cps signal with no distortion. 
Subjecting the thermocouple to various step changes in temperature, it 
was folndothat«for stillliwater , =. 5 see. tor stall air. — ly ecco. 
and for response to the variation of human breath, t = .6 sec. The 


measured time constant of the potentiometer was .9 sec + .2 Sec. 

Based on these observations, a reference temperature of 32.2°F was 
asswned for all measurements and the variation in ice bath temperature al- 
lowed for in the estimated larger possible negative error. 

Assuming a linear temperature transient over a 1 minute period, a vari- 
ation from 190°F to - 50°F would result in a rate change of - 4.0°F/sec. 
Hence, for the non-condensing curves, this effect alone could contribute an 
error of about 5°F, in the thermocouple (t = .6) and potentiometer circuit 
(t = .9 sec). Hence, in the non-condensing curves, an error of + 4.5°F, 

- 1.5°F is estimated. In measurements involving condensation, care was 
taken to let temperature readings reach a relatively constant level before 
data (photos) were taken. No transient effects are estimated for those 


Cases. 





bri), h PATE 4 
- . 45s 


ner 

PRESSURE 

TRMPERARHRE MEASUREMENTS 

Pressure readings taken by the mercury manometer board are estimated 
accurate within + .1" Hg, this would amount to + 1.25% error in the pressure 
ratios near station 18 on the figures for Nozzle I. 

Minimum pressure readings for Nozzle II were in the vicinity of 28" Hg, 
well beyon@ the region of condensation. Hence, pressure ratios in this vi- 
cinity are estimated accurate to within .5%. 

In the region of condensation (station 10) pressure levels range near 
80" He. Hence, error in pressure ratios in this region is estimated at 
+ .125%, two orders of magnitude less than the pressure variation caused by 
condensation at this point. 

Error has been introduced into the values of stagnation pressure wher- 
ever they have not been directly measured by the manometer board owing to 
the high pressure level,or the stagnation tank pressure gauge, as occurred 
several times for Nozzle I. General observations of the behavior of the 
pressure ratio at the throat have shown that although its value remained 
remarkably close to .535 for all pressure levels, it had a slight tendency 
to decrease with increasing pressure. Based on these observations (and 
others with air, where average P/p, measured .523 and the same tendency was 
noted) it is estimated that the pressure ratio at the throat could have been 
lowered to a minimum value of .525 at the highest pressure levels. Hence, 
calculated stagnation pressures are estimated to have an error of + 1%, - 2%. 

Measurements in Nozzle II used the stagnation pressure gauge for all 
determinations of stagnation pressure after the gauge was calibrated and 
found to be accurate to within 1/2 psi at all pressures within the range of 
experimentation. Consequently, all stagnation pressure readings for Nozzle 
If are estimated accurate to within 1/2 psi and static pressure readings 


Wien al ies 





3. 

DENSITY 

Distance between optical glass side walls in the interferometric meas- 
urements was .956 inches. At the wall the boundary layer was about .006" 
at the throat and somewhat larger downstream (See Figure 12). Within the 
boundary layer densities are lower than in the mainstream, due to temperature 
recovery effects. This will cause a small error in the measured interfero- 
metric fringe shift data when applied to determination of density in the mid- 
stream. An upper limit on error in density measurements was estimated by 
assuming a thickness of the order of the displacement thickness (also given 
an upper limit of 6* = 1/3 &) to be of temperature equal to the stagnation 
cemperature. 


Near the throat, assuming 6* = .002" 





26* = ,00}" b= 7956" = Aon 
T 0 
-Sozen — 6 8 ° 6 = e 8 
“§ ° stream 
Hence, 
ap les 20 | also, oe Sees 60 : 


The fringe shift due to difference in density over length & is 


neice 
Hence, 
5 | OF cern Seo ee ae 
S L(p -o_) roe coe 


Similar estimates further downstream indicated that maximum error in 


density measurements probably do not exceed .2%. 














